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Abstract 

The postulational basis of classical thermodynamics has been expanded to in¬ 
corporate equilibrium fluctuations. The main additional elements of the proposed 
thermodynamic theory are the concept of quasi-equilibrium states, a definition of non¬ 
equilibrium entropy, a fundamental equation of state in the entropy representation, and 
a fluctuation postulate describing the probability distribution of macroscopic parame¬ 
ters of an isolated system. Although these elements introduce a statistical component 
that does not exist in classical thermodynamics, the logical structure of the theory 
is different from that of statistical mechanics and represents an expanded version of 
thermodynamics. Based on this theory, we present a regular procedure for calculations 
of equilibrium fluctuations of extensive parameters, intensive parameters and densi¬ 
ties in systems with any number of fluctuating parameters. The proposed fluctuation 
formalism is demonstrated by four applications: (1) derivation of the complete set of 
fluctuation relations for a simple fluid in three different ensembles; (2) fluctuations in 
finite-reservoir systems interpolating between the canonical and micro-canonical en¬ 
sembles; (3) derivation of fluctuation relations for excess properties of grain boundaries 
in binary solid solutions, and (4) derivation of the grain boundary width distribution 
for pre-melted grain boundaries in alloys. The last two applications offer an efficient 
fluctuation-based approach to calculations of interface excess properties and extraction 
of the disjoining potential in pre-melted grain boundaries. Possible future extensions 
of the theory are outlined. 
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1 Introduction 

1.1 Historical background and goal of the work 

Fluctuations of thermodynamic properties play a crucial role in many physical phe¬ 
nomena and diverse applications. Fluctuations are especially important in nanometer-scale 
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systems where they can lead to a large variability of mechanical and functional properties 
and create noises affecting performance of devices. Fluctuations are unavoidable in molecu¬ 
lar dynamics and Monte Carlo simulations of materials, where the system dimensions rarely 
exceed a few nanometers. In fact, in many atomistic calculations, equilibrium properties of 
interest are extracted by analyzing statistical fluctuations of other properties [ij-yj. Exam¬ 
ples include the calculations of elastic coefficients of solid materials from strain and/or stress 
fluctuations in molecular dynamics jH-0] or Monte Carlo [§J, |9j simulations; calculation of 


partial molar properties of solutions from concentration fluctuations [10| 
tion of the interface free energy of solid-liquid [hi, 13 


111 ]: and calcula- 
and solid-sold 14-23| interfaces from 
capillary fluctuations or fluctuations of the interface width. 

Presently, fluctuations are primarily discussed in the statistical-mechanical literature and 
typically in the context of specific (usually, very simple) models. At the same time, the ther¬ 
modynamics community traditionally relies on classical thermodynamics (24-28] which, by 
the macroscopic and equilibrium nature of this discipline, disregards fluctuations and oper¬ 
ates solely in terms of static properties. In fact, the very term “fluctuation” has a temporal 
connotation incompatible with the time-independent character of classical thermodynamics. 

There is, however, a third direction in thermal physics that pursues a generalized form of 
thermodynamics that incorporates fluctuations of thermodynamic parameters around equi¬ 
librium. By contrast to the statistical-mechanical approach, such theories seek to introduce 
fluctuations directly into the thermodynamic framework via additional assumptions, postu¬ 
lates or similar elements of the logical structure. The goal of such theories is to expand the 
scope of thermodynamics by introducing statistical elements while preserving the traditional 
axiomatic approach that distinguishes this discipline. It is this direction in thermal physics 
that constitutes the subject of the present paper [] 

The logical foundations of classical thermodynamics have been the subject of research 
over the past hundred or more years, starting with the works of Caratheodory (29| and 
Ehrenfest (30) and continuing in modern times [2(1 IaI. 3ll 3(J. In spite of the fundamental 
importance of the traditional three laws of thermodynamics, they are essentially a reflection 
of the historical development of the discipline and do not constitute an autonomous and 

27l is de- 


logically complete structure. The four-postulate structure proposed by Callcn 
tached from the historical context, deeply thought-through, but still far from complete. The 
most rigorous and complete postulational basis of classical thermodynamics has been for¬ 
mulated by Tisza [28|. His theory, called the Macroscopic Thermodynamics of Equilibrium 
(MTE), presents an elegant logical structure comprising a set of interconnected definitions, 
postulates and corollaries. This rigor comes at a price: Tisza’s MTE is restricted to a certain 
class of rather simple thermodynamic systems. Nevertheless, it demonstrates an approach 
that can be applied for the construction of similar postulational structures for other classes 
of systems. 

Unfortunately, attempts to create a similarly rigorous thermodynamic formalism that 


Although we do not wish to enter into terminological discussions, we point out that an expanded ther¬ 
modynamic theory that incorporates statistical elements such as fluctuations could be called statistical 
thermodynamics. This term would distinguish it from both statistical mechanics and classical thermody¬ 
namics. Unfortunately, the term statistical thermodynamics is already used with several different meanings, 
most notably as a collective reference to statistical mechanics and thermodynamics (which we refer to as 
thermal physics). 
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would include fluctuations have not been very successful. The thermodynamic fluctuation 


theory by Greene and Callen (371 . 138 ] and Callen’s Postulate II' [26p (which generalizes his 


entropy Postulate II to include fluctuations) turned out to be insufficient. They required 
additional assumptions, such as the approximation of average values of thermodynamic 
properties by their most probable values, and relied on the assumption that there is no dis¬ 
tinction between the “canonical thermodynamics” and “microcanonical thermodynamics”. 
The latter assumption was later criticized by Tisza and Quay (39j from the standpoint of 
statistical mechanics. Tisza and Quay’s [39| own treatment of fluctuations was an amalga¬ 
mation of statistical mechanics and thermodynamics, and in this sense, did not achieve the 
goal either. In fact, the authors j39| admit that their theory, called the Statistical Thermo¬ 
dynamic of Equilibrium, is only a stepping stone toward a future more general theory whose 
development encounters major difficulties. 

There are two major challenges in the development of a thermodynamic theory of fluctu¬ 
ations: (1) appropriate definition of non-equilibrium entropy that would seamlessly connect 
with the system of other definitions and postulates of thermodynamics, and (2) formulation 
of a general fluctuation law of thermodynamic parameters. Historically, the first work ad¬ 
dressing these issues was Einstein’s 1910 paper on critical opalescence j4cj. In his theory, 
the point of departure was Boltzmann’s equation 


S — k In W + const, 


( 1 ) 


where k is Boltzmann’s constant and W is the probability associated with the entropy S (our 
notations are slightly different from Einstein’s). Einstein considered a completely isolated 
thermodynamic system with a fixed energy E. He assumed that all equilibrium and non¬ 
equilibrium (fluctuated) thermodynamic states of the system can be described by a set of 
macroscopic parameters Ai,..., X n . For a discrete set of such parameters, he proposed to 
interpreted W as the fractions of time spent by the fluctuating system in each of the states 
compatible with the given energy E. Based on this interpretation, Einstein inverted Eq. (JT]) 
to obtain the probability of a macroscopic state and thus a given set of parameters Ai,..., X n : 


W = const • exp 


S 


( 2 ) 


For continuos state variables, the probability that the parameters Ai,..., A n be found between 
Ai and Ai + dX\ ,..., X n and X n + dX n is 


dW 


const • exp 


py) 


X * * * ? 


( 3 ) 


where Sq is the maximum entropy. Einstein’s theory does not make any reference to micro¬ 
states of the system. Instead, it relates the entropy directly to the probabilities of macro¬ 
scopic states. As such, this theory belongs more to the realm of thermodynamics than 
statistical mechanics J^] 


2 Callen’s Postulate II' [26;] was abandoned in a later edition of his book [27[ . 

3 The first statistical-mechanical theory of fluctuations was developed by Gibbs [Zlj |. whose work was 
apparently unknown to Einstein at the time. Gibbs’ treatment of fluctuations was based explicitly on classical 
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Equations (J2J) and (J3J) essentially constitute a fluctuation law that can be built into 
the formalism of thermodynamics. The fluctuations described by these equations occur 
in an isolated system, thus the respective entropy has a micro-canonical meaning. For 
canonical systems, the distribution function of macroscopic parameters can be derived by 
considering a relatively small (but macroscopic) subsystem of an isolated system and treating 
the remaining part (complementary system) as a reservoir. 

The situation is complicated by the fact that the most advanced thermodynamic theories 
of fluctuations [26[ 39] postulate the fluctuation law directly in a canonical form and derive 


properties of isolated systems as a liming case of canonical. This approach is different from 
Einstein’s 40] and is rooted in statistical-mechanical arguments. The treatment of a canon¬ 
ical system as an asymptotic case when the complementary system tends to infinity has 
certain conceptual disadvantages in rigorous formulations of statistical mechanics. Alterna¬ 


tive theories [39], [43|, |44] postulate the canonical distribution for an infinite reservoir from 
the outset and derive properties of isolated and finite systems only later. Thus, postulating 
thermodynamic fluctuations in a canonical form is more consistent with the aforementioned 
statistical-mechanical theories. 

In the present paper, we take the view that, if a thermodynamic fluctuation theory is 
to be relatively autonomous, then the simplicity of its internal logical structure is more 
important than consistency with all aspects the formal structure of statistical mechanics. 
We thus propose to build the theory on a fluctuation law postulated for an isolated system. 
The law which will be adopted in this work will look similar to Eq. (|3]h except for certain 
refinements in justification and interpretation. 

We emphasize that this paper is focused on equilibrium fluctuations and will not discuss 
driven or other strongly non-equilibrium systems. Furthermore, given that our approach 
is purely thermodynamic, we are only interested in statistical distributions of fluctua tin g 
parameters without attempting to describe the actual dynamics of the fluctuations j45[ 


1.2 Organization of the paper 

The outline of the paper is a follows. In Secs. [2HZ] we formulate the basic assumptions 
of the theory describing equilibrium fluctuations in an isolated system. We then consider 
the usual quadratic expansion of the entropy around equilibrium and arrive at the Gaussian 
law of fluctuations and a set of relations for mean-square fluctuations and covariances of 
thermodynamic properties (Sec. [8]). In Sec. [9] we derive the fluctuation law for generalized 
canonical systems, treating the system of interest and the reservoir as a combined isolated 
system. We derive the Gaussian law of canonical fluctuations in two forms, called the entropy 
scheme and the energy scheme, and present the fluctuation formalism in both schemes. To 
illustrate the calculation techniques, we compute fluctuations of all properties of a simple 
fluid in three different ensembles. In Sec. m we analyze systems connected to a finite- 
size reservoir. The finite-reservoir ensemble interpolates between the micro-canonical and 
canonical ensembles and has been implemented in recent computer simulations 2jJ, I47H55 


dynamics of particles, with micro-states defined in the momentum-coordinate phase space. Einstein’s theory 
[40( is obviously more general, even though Eqs. (J2|) and @ per se do not suggest any means of calculating 
the probabilities or the entropy for any particular system. For a more detailed comparison of the Gibbs and 
Einstein approaches see Tisza and Quay |39| and a more recent paper by Rudoi and Sukhanov jUJ. 
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This section demonstrates the ease with which the theory can handle fini te-size systems as 
opposed to alternate theories 26, 39} built directly on the canonical formalism. 

Sections EH and EH are devoted to applications of the theory to fluctuations in a single- 
phase interface, called a grain boundary (GB), in a binary solid solution. This is a non-trivial 
problem, which is addressed by conceptually partitioning the system into an imaginary 
perfect crystal (grain) and an excess system representing the GB, each described by their 
own fundamental equation. The goal is to describe fluctuations in the excess system. In 
Sec. EH we develop a set of equations linking fluctuations of some interface properties to 
equilibrium values of other properties, such as the excess heat capacity, excess elastic moduli 
and others. An interesting result here is the prediction of fluctuations of the interface free 
energy, an important excess quantity that is usually treated as a static property. Next, 
we address the problem of width fluctuations of a pre-melted GB represented as a thin 
liquid layer subject to a disjoining pressure. We derive the distribution function of the GB 
width expressed through the disjoining potential describing interactions between the two 
solid-liquid interfaces bounding the liquid layer. While the functional form of the width 
distribution is known for single-component systems (16-18, 20], it has not been previously 
derived for binary mixtures. Finally, in Sec. EH we summarize the paper and outline future 
work. 

Although the paper is focused on thermodynamics, in Secs. El and 19.21 and a few other 
parts of the text we do discuss statistical distributions of micro-states and other aspects of 
statistical mechanics. It should be emphasized that these discussions are not part of the 
proposed fluctuation theory. They are only included here to provide a statistical-mechanical 
interpretation of certain thermodynamic statements or demonstrate their consistency with 
principles of statistical mechanics. In addition, we slightly deviate from the traditional 
macroscopic terminology of classical thermo dynamics and talk about the number of particles 
(instead of moles) and ensembles (instead of types of imposed constraints). Obviously, this 
is only a matter of terminology that does not compromise the thermodynamic nature of the 
theory. 

Finally, we generally prefer a smooth development of the logic and do not follow the 
mathematical style of presentation 


39 in which all statements are broken into 


definitions, postulates, corollaries and other deductive steps. 


2 Definition of equilibrium fluctuations 

Consider an isolated thermodynamic system, i.e., a system incapable of exchanging heat 
or matter with its environment or performing work on the environment. Some physical 
properties of an isolated system are strictly fixed by conservation laws (e.g., the total energy, 
total momentum, total amounts of chemical element^]), whereas other properties can vary. 
Thermodynamics postulates that all non-conserved properties averaged over a certain time 
scale txD-, which we call the thermodynamic time scale, eventually stop varying and remain 
constant for as long as the system remains isolated. The quiescent state of an isolated 
system in which all properties remain constant on the tm time scale is called the state of 
thermodynamic equilibrium. 

4 In a system without nuclear reactions. 
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Even after thermodynamic equilibrium has been reached, non-conserved properties con¬ 
tinue to vary on a shorter time scale, tf C t TD , called the fluctuation time scale. Such 
continual variations of properties of an equilibrium system are called equilibrium fluctua¬ 
tions. 

For any physical property X , let X be its value averaged over the thermodynamic time 
scale tTD- During the equilibrium fluctuations, instantaneous values of X randomly devi¬ 
ate from X. The goal of the fluctuation theory is to predict the probability distribution 
of any property X around its average X and compute the moments and other statistical 
characteristics of this distribution. 


3 The fundamental equation 


Consider an isolated thermodynamic system whose equilibrium properties averaged over 
the tTD time scale are fully defined by a set of n conserved extensive parameters X 1; ..., X n . 
(Tisza |28j refers to such conserved extensive properties as “additive invariants”.) For ex¬ 
ample, all equilibrium properties of a simple fluid enclosed in an isolated rigid box are fully 
defined by its energy E, volume V, and the number of particles nE As long as the pa¬ 
rameters Ai,..., X n remain fixed, the system remains in a state of equilibrium. One can 
temporarily break the isolation, change all or some of the parameters Xi ,..., X n , and then 
isolate and re-equilibrate the system to a new state. By repeating such perturbation/re¬ 
equilibration steps, one can vary the parameters Xi,...,X n and measure all equilibrium 
thermodynamic properties of the system as functions of X \,..., X n . 

In particular, the described procedure can be applied to measure the system entropy S as 
a function of the conserved parameters Xi ,..., X n . Here, the entropy is treated as a property 
defined in statistical mechanics. Namely, for an equilibrium isolated system, S is identified 
with the micro-canonical entropy S = k In Q max , where D max is the maximum number of 
micro-states compatible with the given set of parameters X \,..., X n (see Sec. [7] for a more 
detailed statistical-mechanical interpretation of entropy). The entropy S so defined will be 
referred to as the “equilibrium entropy” in order to distinguish it from the “non-equilibrium 
entropy” defined later. 

The function 

S = S (Xi, ...,X n ) (4) 


is called the fundamental equation of the system in the entropy representation 24, 2tj-28 
More accurately, Eq. (J4J) is the fundamental equation of a particular phase of the system. 
Different phases are specified by different fundamental equations (pfj). For a simple fluid, the 
fundamental equation has the form S = S(E , V, N). 

It is important to recognize that, at this point of the development, the entropy S has 
been defined only for an equilibrium isolated system. Furthermore, the arguments X 1; ..., X n 
of the fundamental equation (J4J) are understood as conserved properties that remain fixed 
(no fluctuations!) in every given state of the isolated system. 


5 This can be considered a definition of the simple fluid. Physically, it is a homogeneous single-component 
fluid without electric, magnetic or other contributions to its thermodynamic properties. 
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4 The non-equilibrium entropy 


When a non-equilibrium isolated system approaches equilibrium, it is often possible to 
define a property called the non-equilibrium entropy. Consider two examples. 

Example 1. When the system is close enough to equilibrium, it can be mentally parti¬ 
tioned into small but macroscopic subsystems that can be considered as isolated and equi¬ 
librium for a period of time on the order of t q . Such subsystems are called quasi-equilibrium 
and the entire isolated system is said to be in a quasi-equilibrium statejj On the quasi- 
equilibrium time scale t q , each subsystem a can be described by a fundamental equation 
S a = S a (Xf, ...,Xf), where ..., Xf are extensive properties assumed to be fixed. The 
non-equilibrium entropy S is defined as the sum of the entropies of all quasi-equilibrium 
subsystems: 

(5) 


s (Ax,..., A m , x u ..., x n ) = J2 S - ■■■> 


Here, Ai,..., A m are so-called internal parameters describing the distribution of the conserved 
extensive properties X \,..., X n of the entire isolated system over its quasi-equilibrium sub- 
systemsQ As the system evolves towards equilibrium, these internal parameters vary on 
the time scale t q , as do the local parameters Xf, ...,X£ of the subsystems and thus the 
non-equilibrium entropy S. Note that the quasi-equilibrium system as a whole does not 
have a fundamental equation; its properties depend not only on the total amounts of the 
extensive quantities Xi, ..., X n but also on their distribution among the subsystems. 

As a specific example, consider an isolated cylinder with rigid side walls, rigid top and 
bottom, and a movable piston dividing the cylinder in two compartments 1 and 2 (Fig. [[]). 
The compartments are filled with different amounts of the same simple fluid obeying a 
fundamental equation S = S(E, V, N ). The system is initially not in equilibrium. Heat and 
particles can diffuse through the piston, but so slowly that each compartment maintains its 
own thermodynamic equilibrium at all times. Thus, as the system drifts toward equilibrium, 
it always remains in quasi-equilibrium. During this quasi-equilibrium process, the total 
energy E = E\ + E 2 , total volume V — V\ + V 2 and the total number of particles N = 
N 1 + N 2 remain fixed. These properties are the additive invariants denoted as X t . The 
non-equilibrium entropy of the system equals 


S{E 1 ,V U N 1 ,E,V,N) = S{E 1 ,V 1 ,N 1 ) + S{E-E 1 ,V-V 1 ,N**N 1 ), 


( 6 ) 


where E\, V\ and N\ play the role of the internal A-parameters describing the distribution 
of the fixed quantities E , V and N between the compartments. 

Example 2. Many systems reach equilibrium with respect to some thermodynamic prop¬ 
erties while other properties are still in the process of equilibration. In such cases, we can talk 
about quasi-equilibrium as equilibrium with respect to some properties without equilibrium 

6 On the time scale t qi the quasi-equilibrium system can be treated as if it was an equilibrium system 
with subsystems separated by isolating walls. 

7 To avoid notational confusion, we note that our parameters Ai,...,A m are independent variables de¬ 
scribing the macroscopic states subject to isolation constraints. In Einstein’s work [4oj, Ai,...,A„ are 
not independent variables. His Eq.([3j) is only valid on the (n — l)-dimensional constant-energy surface 
E( Ai,..., An) = const in the n-dimensional parameter space of Ai,..., \ n . 
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Figure 1: A model of an isolated system: a rigid cylinder divided in two compartments 
containing different amounts of the same substance separated by a movable piston. 


with respect to other properties. For example, chemical reactions usually occur much slower 
in comparison with thermal and mechanical equilibration. There is a certain time scale t q on 
which chemical reactions can be considered as “frozen” and the system can be treated as if it 
were in full equilibrium. Accordingly, the system can be assigned a non-equilibrium entropy 
S computed by ignoring the reactions and treating the amounts of all chemical components 
as fixed parameters. Of course, on the time scale much longer than t q , the entropy slowly 
drifts as the system evolves towards chemical equilibrium. The non-equilibrium entropy so 
defined is a function of not only the conserved parameters Ad, ..., A n , but also some progress 
variables describing the extents of the chemical reactions. The latter thus play 

the role of the internal parameters specifying the quasi-equilibrium states. 

As a simple illustration, consider an isolated box containing a mixture of three chemically 
different gases A, B and C. In the absence of chemical reactions, this mixture is described 
by a fundamental equation S = S(E,V, Na, Nb, Nc). Let the system initially contain the 
amounts of the components N^\ N^ and Nq\ Now suppose that the components can 
transform to each other via a chemical reaction A + B ^ C. Assuming that the reaction is 
slow, the system always maintains mechanical and thermal equilibrium, while its chemical 
composition gradually changes. Let AAA be the change in the number of particles A relative 
to the initial state. The non-equilibrium entropy is defined by 


S(XN a E, V, Ny, AT 


(®) 7 \ 7 "(*) 

B > 



S(E, V, Nf - AN A , N® - AAA, N® + AAA). (7) 


Here, AAA characterizes the progress of the chemical reaction and plays the role of an 
internal parameter A, whereas E,V, N A \N^\ is a set of fixed extensive parameters 
Ah, •••, X n . 


5 The second law of thermodynamics 

In the extended thermodynamics including fluctuations, the second law (entropy postu¬ 
late) can be formulated as two statements: 

1. The non-equilibrium entropy S of an isolated system averaged over the thermodynamic 
time scale tro, which we denote S, increases with time and reaches a maximum value 
S when the system arrives at equilibrium. 














2 . The individual values of S fluctuate but can never exceed the equilibrium value S. 

Note that this formulation of the second law is stronger than the frequently used weak 
formulation. The latter only states that entropy should increase but does not guarantee 
that an equilibrium state will be eventually reached^ Clearly, to talk about equilibrium 
fluctuations we must first ensure that the system actually arrives at equilibrium. 

Once equilibrium has been reached, the average X of any thermodynamic property 
X remains constant while instantaneous values of X can randomly deviate from X on 
the fluctuation time scale tf. This includes fluctuations of the non-equilibrium entropy S. 
However, the latter is a special case because it can only fluctuate down from its maximum 
value S. As a result, the average entropy S is always slightly smaller than S. 

The time evolution of an isolated system approaching equilibrium is illustrated schemat¬ 
ically in Fig. [2j The plot assumes the existence of three different time scaled such that 


txD tf t q . 


( 8 ) 


It should be remembered that S is only defined on the quasi-equilibrium time scale t q . The 
inequality (JHJ) reflects the important assumption of the fluctuation theory that all states 
arising during the fluctuations are quasi-equilibrium states. As such, they can always be 
assigned a non-equilibrium entropy S. The assumption that the approach to equilibrium 
(relaxation) and equilibrium fluctuations realize the same, quasi-equilibrium, states of the 
system bears a certain analogy with the fluctuation-dissipation theorem in statistical me¬ 
chanics 

The non-equilibrium entropy S reaches its maximum value S for some set of internal 
parameters A°, ...,A^ called the equilibrium internal parameters. The latter must satisfy 
the necessary conditions of extremum, 



= 0, i — 1,..., m, 




and the relation 


s (A?,.... A;:„.Y h ...,X n ) = S(X U .... X n ), 


(9) 


( 10 ) 


where the right-hand side is the fundamental equation of the substance. This is illustrated 
schematically in Fig. |3](a) for a single internal parameter A. 

As an example, consider again the isolated cylinder with two compartments discussed 
in Sec. |4] (Fig. [[]). We have three internal parameters: Ai = Ei, X 2 = V\ and A 3 = N\. 
Applying Eq.Q to S given by Eq. (J 6 J) , we have 


dS dS(Ei, Vi, Ni) dS(E 2 ,V 2 ,N 2 ) 


dE 1 


dE 1 


dE 2 


1 1 

“ T\ ~ T2 ~ °’ 


( 11 ) 


8 The weak formulation of the second law of thermodynamics has a “negative” (prohibitive) character 
stating what cannot happen - the entropy cannot decrease; whereas the strong formulation makes a “positive” 
statement predicting what will happen - the system will reach equilibrium and the entropy a maximum. 

9 Generally, the time scales may depend on the observable. For the present discussion, it will suffice to 
represent each time scale by one “characteristic” order of magnitude. 
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S 



time 


Figure 2: Schematic temperature dependence of entropy of an isolated system. The system 
is initially in a highly non-equilibrium state in which the non-equilibrium entropy S is 
undefined. As the system evolves towards equilibrium, it eventually reaches the quasi¬ 
equilibrium stage in which the non-equilibrium entropy S can be defined on a time scale 
t q . The subsequent time dependence of S is shown by a broken line to emphasize that 
individual values of S can only be measured within time intervals on the order of t q . The 
values of S can never exceed the maximum (micro-canonical) entropy S. By the second 
law of thermodynamics, S averaged over the thermodynamic time scale txD (S, shown by 
the dashed blue line) monotonically increases with time and eventually levels out when the 
system reaches equilibrium. In the equilibrium state, S can randomly deviate down from 
its maximum value S on the fluctuation time scale tf such that t q <C tf <C txD, always 
returning to the maximum value when the fluctuation is over. 




Figure 3: Schematic diagram of equilibrium fluctuations in an isolated system with one 
internal parameter A. (a) Non-equilibrium entropy S reaches its maximum S at the equilib¬ 
rium value of the internal parameter A 0 . The arrows symbolize equilibrium fluctuations of A 
on either side of A 0 accompanied by downward deviations of S from S. (b) The probability 
density of A has a sharp maximum at A = A 0 of a height W m . 
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( 12 ) 


as as(Ei, Vi, Ni) 


dVi 


dVi 


dS(E 2 ,V 2 ,N 2 ) _ Pl p 2 _ 

dV 2 T, T 2 


dS(E 1 , Vi, Ni) dS(E 2 ,V 2 ,N 2 ) 


-)- — — o 
T, T 2 


(13) 


aw aw dw 2 

where we used the standard thermodynamic relations (2514271. f58| dS/dE = 1/T, dS/dV = 
p/T and dS/dN = —/i/T (T being temperature, p pressure and /i chemical potential). We 
thus recover the well-known thermodynamic equilibrium conditions T\ — T 2 , p\ = p 2 and 
Pi = p 2 . These conditions must be satished for the equilibrium parameter set (E®, V®, N [). 

As another example, the non-equilibrium entropy of a gas mixture with the chemical 
reaction A + B ^ C is defined by Eq.(I71) with a single internal parameter A = AN a- The 
necessary condition of equilibrium reads 


dS(E, V , W« - AN a, N%> - AN A , N$> + AN a ) 

8AN a 

dS(E, V, N a , N b , N c ) dS(E, V, N A , N B ,N C ) dS(E , V, N A , N B , N c ) 


0 


rW 


3iV A 


<9W 


B 


aw 


c 


= 0. 


(14) 


We obtained the standard chemical equilibrium condition p A + Pb = Pc, which must be 
satisfied at some A 0 = AIV^. The equilibrium amounts of the chemical components are then 


K = n a 


(*) 


AiV^, W° = wf } 


B 


AN^ and W° = N$ 


A N° a . 


6 The general law of fluctuations 

During equilibrium fluctuations, the internal parameters randomly deviate from the equi¬ 
librium values Xi ,..., A° n on the time scale oft/. These fluctuations of the internal parameters 
are accompanied by deviations of the non-equilibrium entropy S down from its maximum 
value S [Fig. [3](a)] . 

The Fundamental Postulate of the thermodynamic theory of fluctuations states that the 
probability density W(Ai,..., A m ) of internal parameters of an equilibrium isolated system is 


W(Ai,...,A m ) = A exp 


‘S'(Ai,..., A m ) 

k 


(15) 


where A is a constant. To simplify the notations, we do not display the fixed parameters 
Xi,X n as arguments of W and S. 

By definition, W (Ai,..., X m )dXi...dX m is the probability of finding the internal parameters 
in the infinitesimal region dX\...dX m of the parameter space. The constant A can be found 
from the normalization condition 


11 (A^,..., A m )dAj.. ,dX m 1, 


(16) 
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where the integration extends over the entire parameter space. Since S reaches its maximum 
value 

S = S( a;,...X) (17) 

at the equilibrium parameter set A?,..., A^ [cf. Eq. ffTOjl ]. W has a peak of the height 


W m 



(18) 


at A°, A^ [Fig. (3](b)]. In other words, A°, •••, A^ is the most probable parameter set. Using 
this fact, Eq. (TTTTI) can be recast as 


VF(A 1 ,..., A m ) 


W m exp 


S{\i, •••, A m ) — S 
k 


(19) 


Equation (1T9|) is the central equation of the thermodynamic fluctuation theory. It is 
the departure point for thermodynamic analysis of all statistical properties of equilibrium 
fluctuations. 

The average (expectation) value of any internal parameter A, is obtained by 


A,: = W„ 


A i exp 


^(Ai,..., A m ) - S 
k 


d /\ ^... 


( 20 ) 


Note that this average does not generally coincide with the most probable value A°. It 
is common to use the symbol A to denote deviations of fluctuating properties from their 
average values. For example, 

AAi = Xi- A*. (21) 


The covariance (second correlation moment) of any pair of internal parameters can be com¬ 
puted by 


AXiAXj 


W m 


I A \ A \ ovn 

'S(Xi, ..., A m ) - s' 

I l -A / \ r^L - \/\j LAjJ 

l 

_i 


dX\. ..dX m , 


( 22 ) 


while the mean-square fluctuation of A* is 


(AA,)^ 


W m 


/ f A A A 2 pYn 

'^(Ai,..., A m ) - s' 


l 

_i 


dXi...dX r 


(23) 


7 Statistical-mechanical interpretation of fluctuations 

In this section we will discuss the general law of fluctuations formulated above from the 
viewpoint of statistical-mechanics. 

As before, we consider an isolated system whose equilibrium state is defined by a set of 
fixed parameters X \,..., X n . If multiple observations are made, the system can be found in 
different macro-states, each specified by a set of internal parameters Ai,..., X m . Let us first 
consider a system with discrete A-parameters. Each parameter set Ai,..., A m is implemented 
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by a certain number O of micro-states. The latter number is sometimes called the degeneracy 
of the macro-state. The macro-states are analogs of the equilibrium and quasi-equilibrium 
states introduced in the thermodynamic theory (Sec. [4]). The degeneracy of a given macro¬ 
state is a function of the internal parameters Ai,..., X m and the fixed parameters X \,..., X n : 

0 = 0 (Ai,..., A m , Xi, ...,X n ). (24) 


It is postulated that for an isolated equilibrium system, there is one unique macro-state 
for which the degeneracy reaches a maximum value Ft max . Denoting the internal parameters 
corresponding to this macro-state A°, ...,A^, the maximum degeneracy depends solely on 
the fixed parameters X 1} X n : 


(AV .... X„) = SI (a;, .... A;:.. A-,..... A'„) . 


(25) 


For a given set of fixed parameters X t , ..., X n , the degeneracy O of an equilibrium isolated 
system fluctuates with time from its maximum value 0 ma3 . down to smaller values. These 
equilibrium fluctuations follow the Fundamental Postulate of Statistical Mechanics stating 
that the probability P of a single micro-state is the same for all micro-states and depends 
only on the fixed parameters X \,..., X n : 


P = P(X U -, X„). 


(26) 


This statement is often referred to as the micro-canonical distribution. It immediately 
follows that the probability of any macro-state is simply PD. Thus, macro-states with a 
larger degeneracy have a proportionately higher probability to occur. The macro-state with 
the largest degeneracy fl max is the most probable one. 

Analysis of fluctuations relies solely on the fact that P is the same for all micro-states 
and does not require a specific knowledge of P. We note, however, that P can be found 
from the normalization condition stating that the sum of probabilities of all possible macro¬ 
states of the system be unity. In the approximation where fluctuations are neglected, only 
the most probable macro-state with the highest degeneracy Q max is considered. Then the 
normalization condition reads PFt max = 1, giving 


P 


1 


O 


max 


(27) 


On the other hand, when fluctuations are significant, all possible degeneracies from 0 = 1 
to Pl max must be included. From the normalization condition 


o 


max 


E™ = 1 

0=1 


(28) 


we obtain 


P = 


Xlma X (Xlma X T 1) 


which for large systems can be approximated by 2 /Xlf nax 


( 29 ) 
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Regardless of the specific value of P, the product Phi is the probability of a given macro¬ 
state and thus the probability 


W Xt .A„ = ™< Ai,...,A m ) 


(30) 


of the corresponding set of internal parameters. Here, as in Sec. [6] we simplify the notations 
by suppressing the fixed parameters X n as arguments of P and fl The probabil¬ 

ity W \reaches a maximum value, denoted W m , at the most probable parameter set 
A°,..., A°j, for which 

W m = PQ (A?,..., X° m ) = PVL max . (31) 

Combining Eqs. (l30]i and f|3Tjl . 


W; 




= W r . 


(Ai,..., A m ) 
1 H 

- - rn n:r 


(32) 


We now define the non-equilibrium entropy of an isolated system by 


S (Ai,... ? A m ) — k In (Ai,A m ) 


(33) 


and the equilibrium entropy by 

S = k In O max . (34) 

It immediately follows that the non-equilibrium entropy can never exceed its equilibrium 
value S. It also follows that in equilibrium, S can occasionally fluctuate down from S but 
eventually returns to S as the most probable value. 

The probability of any given set of internal parameters given by Eq. ((32]) can now be 
rewritten in the form 


H / A 1 ,...,A m = W m ex p 


‘S'(A,,..., A m ) - S 


(35) 


This equation treats Ai,...,A m as discrete parameters. For continuous parameters, Eq. (j35j) 
remains valid but gives the probability density of the internal parameters. The obtained 
Eq. (j35]l is identical to Eq. (ll9]l of the thermodynamic theory of fluctuations (Sec. [6]). 


8 The Gaussian law of fluctuations 


In macroscopic systems, the peak of W at the equilibrium parameter set A°, ...,A^ is 
extremely sharp. In most cases, one can safely replace the exponent in Eq. fllhl) by its Taylor 
expansion at A?,..., A^ truncated at quadratic terms: 


i=l 




A?) (A, 


A?), (36) 


where the subscript 0 indicates that the derivatives are taken at A?, •••, A^. The linear terms 
vanish by the extremum condition ([9]). As a result, Eq. (fl9|) reduces to the multi-variable 
Gaussian distribution 


W(Xi ,..., A m ) = W rn exp 




i,j=1 


(37) 
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with the symmetrical matrix 


A h - k 


d 2 S 

d\id\j 


(38) 


which we call the stability matrix of the system. For normally stable systems this ma¬ 
trix must be positive-definite0 Applying the normalization condition (fT 6 |h the Gaussian 
distribution can be rewritten as 





m 

~2 A b(A * 

i,j =1 


A“)(Aj 



(39) 


where A is the determinant of A,-,- 

We can now make use of the well-known mathematical properties of the normal distri¬ 
bution. The average value of any parameter A$ coincides with its most probable value, 


Aj = A°, i = l,...,m. 


(40) 


For covariances of internal parameters we have 

AXiAXj = A"- 1 , i , j = 1, ...,m, (41) 

where Ah 1 is the inverse of the stability matrix A ir In particular, the mean-square fluctua¬ 
tion of a parameter A i equals 


(AA i ) 2 = A ii 1 , i — 1,..., m. (42) 

We can also introduce the following linear combinations of internal parameters: 

m 

<Pi = y^AjjAj, i — 1,..., m. (43) 

j =i 

It can be shown that 

Aipi = 0, i — 1 ,..., m, (44) 

AifiAXj = Sij, i,j = 1,..., m, (45) 

AipiAipj = A i:j , i,j = (46) 

and therefore 

(A^) 2 = An, i — 1, ..., m. (47) 

Thus A ij is the covariance matrix of the parameters (p^. 

An important property of Eqs. flTTf . (1421) . (1461) and (1471) is that their left-hand sides 
contain fluctuations, whereas the right-hand sides are computed in the state of equilibrium 

10 We follow the classification [28| dividing all systems into normally stable, critically stable, unstable and 
metastable. In normally stable systems, all fluctuations are finite. 

11 The Gaussian distribution (IdUl) can be interpreted as a consequence of the central limit theorem for 
a system composed of many statistically independent subsystems, see for example Chapter V of Khinchin 
441. 
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and can be obtained by the standard formalism of equilibrium thermodynamics. In other 
words, these equations relate fluctuations to equilibrium properties of the system. 

As a simple illustration, let us revisit the isolated system divided in two compartments 
(Fig. [[]). Suppose the system is in equilibrium. Consider a particular type of fluctuations 
in which the compartments can only exchange energy at fixed volumes and numbers of 
particles. Thus, the only internal parameter is Ai = E x . The stability matrix (l38|) reduces 
to the scalar 


A 


1 f d 2 S\ 

1 WW 


(48) 


From Eq. flTTTh 


dS _ 1 1 

~dE\ ~ T\~ T~2 


(49) 


and thus 


A = 


<9T\ 1 

+ 


dTo 


k \T 2 dE 1 T 2 dE 2 


kT^Cy 


1 1 
iVi + T 2 


where Tq is the equilibrium temperature of both compartments and 


(50) 


Cv -N\gT) v „ 


(51) 


is the specific heat of the substance per particle at the equilibrium temperature T 0 . Applying 
Eq. (l42|i we obtain the energy fluctuation of compartment 1: 



Vi.iVi 


1 

A 


kT^ c v 


N\N 2 
N\ T N 2 ' 


(52) 


When the second compartment is much larger than the first, N 2 As> Aq, this equation 
reduces to 

f(AE^) = iW 0 2 cv (53) 

\ /Vi ,Ni 

In this limit, the second compartment serves as a heat bath (reservoir) for the first. The 
relative fluctuation (variance) of E\ becomes 



FoM ) 1 / 2 


(54) 


where £ = Ei/Ni is the equilibrium energy per particle. This relation shows that the relative 
fluctuation of energy decreases with the system size as 1/ y/N. 


9 Fluctuations in canonical systems 

9.1 General relations 

We now return to the general fluctuation law formulated in Secs. [6] and will apply it to a 
relatively small (but still macroscopic) subsystem of an equilibrium isolated system. We will 
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ISOLATED SYSTEM 



Figure 4: A subsystem chosen inside an isolated system characterized by a set extensive 
properties X 1; X n . These properties are partitioned between the subsystem and the com¬ 
plementary system (reservoir). The partitioning is specified by the parameter sets X 1 ,X n 
and X[,... ,X ”, respectively. 

adopt the terminology in which we refer to the selected subsystem as the “canonical system”, 
or simply “system”, and to the remaining part of the isolated system as the “complementary 
system” or “reservoir”. 

An example of the situation was already given by the system with two compartments 
(Fig. [T]) when the second compartment is much larger than the first. We will now con¬ 
sider a more general case in which the entire isolated system is characterized by a set of 
conserved extensive properties (additive invariants) Xi,...,X n . These properties are parti¬ 
tioned between our system, Xi,...,X n , and the reservoir, X[,..., X r n (Fig. [4]). For some of 
these properties, the partitioning is fixed once and for all, while m remaining properties can 
“flow” between the system and the reservoir as a result of equilibrium fluctuations. Such 
fluctuations are subject to the conservation constraints 

Xi + X[ — Xi — const, i = 1,..., m. (55) 

We assume that the exchange of the quantities X* occurs slowly enough that we can 
consider both the system and the reservoir as quasi-equilibrium systems. Accordingly, they 
obey the fundamental equations 

S — S (X 1; ..., X n ) 

for the system and 

S r = S r (X[,...,X*) (57) 

for the reservoir. Note that these two equations are generally different, allowing the system 
and the reservoir to be composed of different substances. 

In the subsequent analysis, a special role will be played by the entropy derivatives 

f . = 9S_ 

dXi' 

called the thermodynamic “forces” conjugate to the extensive variables X,. Similar deriva¬ 
tives dS r /dXl define the forces F[ in the reservoir. For example, for a simple fluid with the 


(58) 
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fundamental equation S(E,V, N), the thermodynamic forces are F 1 = 1/T, F 2 = p/T and 
F 3 = —p/T. Note that the thermodynamic forces are intensive variables. 

We define the non-equilibrium entropy of the entire isolated system as 

S(X ll ...,X mi X m+1: ...,X ni X 1 ,...,X n ) = S(X 1 ,...,X n ) + S r (^X 1 -X 1 ,...,X n -X n ^ , (59) 

A 

where X l} ..., X m play the role of the fluctuating internal parameters Aj, whereas the remain¬ 
ing parameters X m+1 , X n , Xj,..., X n remain hxed. The equilibrium value of this entropy 
is 


£ = s(x°,...,x°,x m+1 ,..,x n ) 

+ s r (x, - x°, ...,x m - x° ml x m+1 - X m+1 , ...,X n - X n ) , (60) 

where X°,..,X° are the equilibrium values of the fluctuating parameters. To simplify 
the notations, from now on we will suppress the hxed parameters in all equations, writing 
Eqs. (T59|) and (1601) in the form 


X m XfjiJ , 

( 61 ) 


( 62 ) 


Before analyzing fluctuations, we will apply the equilibrium conditions (j9j) to the internal 
parameters A i = Xi (i — 1, ...,m): 


dS 

dX 


dS 


dX, 


*/ o 


dS r 


dX] 


= Ff-(F[)° = 0, i = l ,...,m, 


( 63 ) 


i / 0 


where, as usual, the symbol 0 indicates that the derivatives are computed for the equilibrium 
parameters A") 1 , ...,X/ n . We thus conclude that in equilibrium, the system and the reservoir 
have equal values of the thermodynamic forces conjugate to the fluctuating parameters. 

We now turn to fluctuations. The probability density of the internal parameters 
Xi, ...,X m is given by Eq. (fT9lh which takes the form 


W(X u ...,X m ) = W m exp 


= W m exp 


S(X 1 ,...,X m )-S 

k 


S(X ll ...,X m )-S(Xl...,X/ n ) 

k 


x exp 


VfV - X u .... X m - X m ) - SAX, - X«, .... X m - X»J 


( 64 ) 


In the second exponential factor, we will expand S r around the equilibrium parameters 
wo yo. 


Sr(X\ - x,,..., x m -,x m ) - 5V(X, - X",..., X m - X°) = - ^/’"(W-X' 1 


2=1 


E 

»j=i 


d 2 S r 


dX/DX], 0 


(X i -X“)(X j -X.9). 


(65) 
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As will be shown later, the second derivatives appearing in the quadratic terms scale as 
1/N r , where N r is the number of particles in the reservoir. Since the reservoir is assumed 
to be much larger than our system, the quadratic terms in Eq. (I65|) can be neglected in 
comparison with the linear terms. Neglecting the quadratic terms constitutes what can be 
called the “reservoir approximation”, in which 

m 

Sr(Xi - x u ..., - X m ) - SriXi - X°, ..., X m - X°J = F?{Xi - X?). (66) 

i =1 


Inserting this equation in Eq. fj64l) . we obtain 

m 

S(X U X m ) - S( XIXI) - V F«(X, - Xf ) 


W(X h ...,X m ) = W m exp 


2—1 


(67) 


Note that this equation contains only properties of the canonical system and does not depend 
on properties of the reservoir. The role of the reservoir is only to impose the thermodynamic 
forces Ff. 

Equation (167|) is an important relation that will serve as the starting point for deriving 
all fluctuation properties of canonical systems. 

In some canonical systems, there can be additional internal parameters, say Yi, Y 2 ,..., Yj, 
that can fluctuate without affecting the reservoir parameters X[,..., X”. In such cases, the 
set of internal parameters is Ah, Y x , Y 2 , ...,1) and Eq. (l6Tf is replaced by 


S(X 1 ,...,X m ,Y 1 ,Y 2 ,...,Y l ) = S(X 1 ,...,X m ,Y 1 ,Y 2 ,...,Y l ) + S r [X 1 -X 1 ,...,X m -X m j . 

( 68 ) 

The equilibrium conditions with respect to the parameters decoupled from the reservoir is 


ai] = (ds_\ =F o =n 

QyA [dYj Q - ’ 


(69) 


where i = 1,2,...,/. It is easy to see that the additional parameters can be incorporated into 
the parameter list X Xl ..., X m by simply assigning them zero thermodynamic forces. Accord¬ 
ingly, Eq. (16711 remains valid, except that the sum with F , ; ° excludes the terms corresponding 
to the decoupled parameters Y 1 , Y 2 ,..., Y{. 


9.2 The generalized canonical distribution 

Before proceeding with applications of Eq. (l6T]) . we return for a moment to the statistical- 
mechanical interpretation of fluctuations discussed in Sec. [3 For discrete internal parameters 
Ah,...,Al m , Eq. (167|) gives the probability Wx 1 ,...,x m °f a given parameter set X 1 ,...,X m . 
Recall that on the quasi-equilibrium time scale t q , our system can be treated as if it were 
equilibrium and isolated with fixed values of X Xl ..., X rn . Thus, its fundamental equation 
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S(Xi ,..., X m ) gives the equilibrium entropy S = k In Q max . The probability P of a single 
micro-state can be obtained from the relation Wx 1 .....x m = which gives 


P 


W 


klrnnT 


W Xl ,...,x m ex P 


k 


Combining this relation with Eq. (l67l) we obtain 


P 


W rn exp 


s(x?,...,x^) + '£,f?(Xi 

i=l 


k 



(70) 


( 71 ) 


which can be rewritten as 


P = A exp 



where 


A = W m exp 


i= 1 


k 


(72) 


(73) 


is a constant. 

Equation (172]) constitutes the generalized canonical distribution describing several sta¬ 
tistical ensembles!^! 

Assuming that one of the extensive parameters, say Xi, is energy, Eq. d72l) can be rewrit¬ 
ten in the form 

/ m \ 

E + ^Toi^O 

i =2 


P = A exp 


kTn 


(74) 


As an illustration, suppose our system is a simple fluid. Either one or two of its extensive 
parameters E , V and N can fluctuate with the third parameter fixed (otherwise the system 
is undefined). For example, suppose only energy can fluctuate while V and N remain fixed. 
From Eq. (l74|) . the probability of a micro-state is 


P = Aex P (—Jr), fixed V, N , (75) 

12 The way we arrived at Eg. (1721) is not intended to be a rigorous statistical-mechanical derivation of 
the canonical distribution. The goal was only to show that the thermodynamic fluctuation relation (1671) is 
consistent with the distribution of micro-states known from statistical mechanics. 
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a relation which is known as the canonical, or NVT, distribution. If both energy and volume 
are allowed to fluctuate, we have X 2 = V with F 2 = p/T, giving the NpT distribution 


P 


A exp 


( E + Po V \ 

V kT, )' 


fixed N. 


(76) 


If the energy and the number of particles can fluctuate at a fixed volume (open system), 
we have X 2 — N and F 2 = —p/T. Equation (l74l) gives the grand-canonical, or pVT , 
distribution 

P — A exp (k — , fixed V. (77) 

Of course, the pre-exponential factor A is different in all three distributions and is related 
to the respective partition functions through the probability normalization condition. 

Note that the distributions (l75lh (|76l) and (1771) contain the temperature T 0 , pressure 
Po and the chemical potential p o corresponding to the exact equilibrium with the reservoir. 
These properties are fixed parameters imposed by the reservoir, by contrast to the fluctuating 
temperature, pressure and chemical potential inside the system. 

We re-emphasize the important difference between Eq. (|72]) and the previously derived 
Eq. (l67l) : Eq. (l72 ] ) gives the probability of a single micro-state of the system, whereas Eq. (l67|) 
is the probability distribution of the parameters X l5 ..., X rn characterizing different macro¬ 
states. 


9.3 The Gaussian law of canonical fluctuations 

We now return to the probability distribution of the fluctuating parameters given by 
Eq. (j67l) . Let us expand the entropy of the system around its equilibrium value, keeping only 
linear and quadratic terms: 


S(Xu .... X m ) - S(X°, ...,X^) = J2 F,°(Xi 

i= 1 



d 2 S 

dXidXj 


(Xt-XfKXj 



(78) 


Inserting this expansion in Eq. (167)) . we arrive at the multi-variable Gaussian distribution 


W(X 1 ,...,X m ) = W m ex p 


i 




(X, - X°)(Xj 


V») 


(79) 


with the stability matrix 


Aij 


1 / d 2 S \ 1 (dF t \ 

k [dX.dXjJo k \dxj 0 


(80) 


In this approximation, the average values X * are numerically equal to the most probable 
values X/. Using the standard properties of the Gaussian distribution (Sec. [8]), we obtain 
the fluctuation relations 


AXiAXj 


A 


-i 

ij > 


hi = 1 


(81) 
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and thus 

(AX,) 2 = A^ 1 , i = 1, m. 
We can also introduce the parameters 

m 

<Pi = jjXj, i = 

3 = 1 


(82) 


(83) 


for which 


<Pi-<Pi = 


1 m 


3 =1 


m 

dX, 




x 3 ) = - k : 


(84) 


Here, AT) is the deviation of F % from its equilibrium value evaluated in the quadratic ap¬ 
proximation (l78lh Using the properties ()46li and ()47j) of the Gaussian distribution, we obtain 
the fluctuation relations for thermodynamic forces: 


AXiAFj = —kdij, 
AFiAFj = k 2 Aij : 

and therefore 

(A W = k 2 An, 


i,j = 1,..., m, 

(85) 

i,j = 1 ,-,m, 

(86) 

i = 1,..., m. 

(87) 


Equations (|82|) and (1871) show that for normally stable states, the derivatives ( dXi/dFi) 0 
and (dFi/dXi) 0 must be positive. At critical points one of these derivatives tends to infinity 
and the respective fluctuation diverges (unless the reservoir has a finite size, see Sec. ITOl) . 

Another form of the foregoing equations is obtained by introducing the thermodynamic 
potential 

m 


e(Ei ,...,F m ) = S-J2 F iXi 


as the Legendre transform of the entropy (Massieu 


i =1 


function) [25- 27|. Then, 


( 88 ) 




<90 

m 


(89) 


and 



/ < 9 2 0 \ 

\WdFjJo 


(90) 


For normally stable states (d 2 Q/dF 2 ) 0 > 0 for all i = 1, ...,m. 

It is interesting to evaluate canonical fluctuations of entropy relative to its equilibrium 
value S(X®, ...,X^). Using the quadratic approximation (178)1 . we have 


m m 

AS = Y, F i Ax i A -,A,Y,A.V r 

i=l i,j =1 


(91) 


Averaging this equation over the distribution ([79]) and taking into account Eq. flSTf . we 
obtain 


m 1 m 

AS = ^TifAW - -k Y AijAXiAXj 

i =1 hj = 1 



(92) 
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Thus, in the Gaussian approximation the canonical entropy equals 

S = S(X°,...,X°)-±k. (93) 

Considering that entropy scales in proportion to the total of number of particles N, the 
relative fluctuation A S/S oc N _1 . Using the same Gaussian approximation, the mean- 
square fluctuation of entropy scales in proportion to N, 

m m 

(KSY = F ° F f AJf.AX, = P°F °AV oc N. (94) 

i,j=1 i,j=1 

where we neglected higher-order moments of the distribution. Thus, the relative fluctuation 
, _\ 1/2 

((AS 1 ) 2 ) /S oc N~ 1/2 . Both estimates show that, with increasing size of the system, 

the canonical entropy S tends to the micro-canonical entropy S satisfying the fundamental 
equation ([4|). While the two entropies are conceptually different, they converge to the same 
numerical value in the thermodynamic limit. 


9.4 Canonical fluctuations in the energy scheme 

The treatment of fluctuations based on the fundamental equation (|4|) is called the entropy 
scheme. In practical situations, energy is always one of the extensive parameters fluctuating 
between the canonical system and the reservoirSuppose X 1 = E and thus F x = 1/T. 
Equation (|4|) can be inverted to obtain 

E — E (Z x ,..., Z n ), (95) 


where Z x = S and Z % = A,, i = 2,..., n. This equation is called the fundamental equation in 


the energy representation, or the energy scheme 28]. In the energy scheme, the thermody¬ 


namic forces Pi = dEfdZi conjugate to the extensive arguments Z x ,..., Z n are P x = T and 
Pi = — TFi, i = 2,..., n. For example, for a simple fluid E = E(S, V, N ) and thus P x = T, 
P 2 = — V an d P3 = /i. 

In many applications, we are more interested in fluctuations of the P-forces rather than 
P-forces. The formalism presented below expresses canonical fluctuations in terms of the 
variable sets Z x ,Z n and P x , ...,P n corresponding to the energy scheme. 

As before, we consider a system coupled to a reservoir and described by n extensive 
properties X x , ...,X n . Out of them, m properties X x ,..., X m can flow between the system 
and the reservoir while the remaining properties X m+ \,..., X n are fixed. Assuming X x = E , 
we transform the variable set X ±,..., X m to Z Xl ..., Z m and perform a Legendre transformation 
to obtain the thermodynamic potential appropriate for the the energy scheme: 


m 

$(Pi,..,P m ) = E-J2 p iZi- (96) 

2=1 

13 While “adiabatic ensembles” in which the system is adiabatically isolated but can still exchange particles 
with the reservoir are helpful theoretical constructions, their practical implementation in experiments or 
atomistic simulations is unfeasible. 
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This potential has the property 


< 94 > 

m 


= -Zi. 


(97) 


Returning to the parameter distribution (j6T|) . we can now recast it in the form 

m 

E-E 0 -Y,lf(Z,-Zf 


W = W m exp 


i —1 


kTr, 


= W m exp 


n 


kTr, 


(98) 


where we denote 


K^E-Eo-Y&iZi-ZO). 


i=l 


(99) 


It can be shown that 1Z has the meaning of the reversible work required for the creation 
of the fluctuated state of the canonical system at a fixed value of the total entropy of the 
system and the reservoir 24], [58 . 


Equation (1981) is the canonical distribution function of the extensive parameters 
Z i, Z m . For small fluctuations, we can expand E — E 0 in these parameters and neglect 
all terms beyond quadratic to obtain the multi-variable Gaussian distribution 


W(Z U ..., Z m ) = W rn exp 
where the matrix 


1 


£ A« (z, - Z«)(Z, - Z») 


*,j= 1 


1 kT„ \dz,az 


3/0 


kTr, \dz 


dP, 


can be called [28] the generalized stiffness matrix. Its inverse 

' dZi\ a 2 <h 


K~ x = kT 0 


dP, 


= ~kT 0 
o dPidPj 


( 100 ) 


( 101 ) 


( 102 ) 


can be then called 28| the generalized compliance matrix. 

Using the standard properties of the Gaussian distribution (Sec. [ 8 ]) , we obtain at once 
the covariances and mean-square fluctuations of all extensive and intensive properties of the 
canonical system in the energy scheme: 


AZiAZj = K tJ 1 , i,j = 1,..., m, 

(103) 

(AZ l ) 2 = A4f i , i — 1,..., m, 

(104) 

AZiAPj = kT 0 5ij, i,j — 1 ,..., m, 

(105) 

APiAPj = ( kToYKij , i,j = 1 ,..., m, 

(106) 

(APi) 2 = (kT o y 2 Ku, i = 1 ,..., m. 

(107) 
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Having these fluctuation relations, it is straightforward to find the covariance of any 
two thermodynamic properties y = y(Zi ,..., Z m ) and z = z(Z 1 ,Z m ). Indeed, for small 
fluctuations, 



(108) 



(109) 


from which 



Similarly, for any two functions y = y(P \,..., P m ) and z = z(P \,..., P m ), 



The same calculation method can be applied when y is a function Zj, ..., Z. m while z is a 
function of P 1 ,..., P m , or when both are functions of mixed sets of Z’s and P’s. 

9.5 Canonical fluctuations in a simple fluid 

To demonstrate the foregoing formalism for computing canonical fluctuations, we will 
apply it to a simple fluid coupled to a reservoir. We will consider three different ensembles, 
deriving equations for the mean-square fluctuations and covariances of all thermodynamic 
properties in each case. The calculations will be conducted in the energy scheme, although 
the entropy scheme could be applied just as well. We are working with the extensive and 
intensive parameter sets (Zi, Z 2 , Z 3 ) = ( S,V,N ) and (Pi,P 2 ,P 3 ) = (T, — p, p), respectively. 
To simplify the notations, we will suppress the index 0 but it is implied that all derivatives 
are taken at the equilibrium state. 

While a number of fluctuation relations for fluids can be found in the literature, in this 
paper we present a complete set of such relations for each ensemble. For the NVT ensem¬ 
ble, there are 5 fluctuating parameters (E, S,T,p, p). Out of 25 covariances that can be 
formed among them, only 15 are distinct due to the symmetry of the covariance. Similarly, 
in the NpT and pVT ensembles, there are 6 fluctuating parameters (E, S,V,T,p, p) and 
(E, S, N,T,p, p), respectively!^! Taking into account the symmetry of covariances, 21 dis¬ 
tinct fluctuation relations exist for each of these ensembles. All these fluctuation relations 
will be derived below in a systematic manner enabled by the proposed formalism. 

The equilibrium properties appearing in right-hand sides of the fluctuation relations will 
be presented in simplest possible form. They can be further transformed to many other 

14 Recall that all intensive properties of a canonical system are allowed to fluctuate, including those which 
are imposed by the reservoir. For example, in the yVT ensemble, the reservoir imposes a certain temperature 
To and chemical potential yo\ however, the actual temperature and chemical potential inside the system 
fluctuate around the imposed values. 
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forms expressing them through specific heats, moduli or other experimentally accessible 
properties. We do not pursue such transformations as they he outside the fluctuation topic 
and depend on the problem at hand. 


9.5.1 NVT ensemble 


Suppose the volume and number of particles in the system are fixed while the energy and 
entropy can fluctuate (canonical ensemble). The stiffness and compliance matrices reduce 
to the scalars 


1 /<9T\ 1 

~kT 0 \ds) VN ~Nte v 


( 112 ) 


and 


K[ i 1 = Nkc v . 


(113) 


By Eq.tHED, 

(AS) 2 = K[ i 1 = Nkc v . 


(114) 


Applying Eq. flllOl) . we find the energy fluctuation 


(A E) 2 



ATi 1 


V,N 


NkTyC v 


and the energy-entropy covariance 


AEAS 



K 


-l 

li 


V,N 


NkT 0 c v , 


(115) 


(116) 


where we used the relation (dE/dS)y,N = T 0 . For temperature fluctuations, we follow 
Eq. (11071) to obtain 

(AT) 2 = (kTofKn = (117) 

We next derive fluctuations of the chemical potential and pressure. Both are dependent 
variables that can be expressed as functions of temperature and the number density of 
particles, p = N/V. The latter is fixed, leaving only T as an independent argument. 
Applying Eq. dllip we obtain 


(Ap) 2 = (kT 0 ) 2 


(A p) 2 = (kT 0 ) 2 
AM? = (kT 0 f (jE) 
ApAT = (kT 0 ) 2 



kT[ 

Nc v 

kT[ 

Nc v 

kT[ 

Ncv 

kTl 

Nc v 



( \ ( 
\&r) p \&r) 



P 


(118) 

(119) 

( 120 ) 
( 121 ) 
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ApAT =(^(§)^ 1 = ^(§) ; ( 122 ) 

For “cross-fluctuations” between extensive and intensive properties, we first apply 
Eq. (11051) to obtain 

AS AT = kT 0 . (123) 


The remaining “cross-fluctuations” are easily computed by the method outlined by Eqs. (jl08j) - 

innj: 

AEAT = (A) AS AT = kT, J, (124) 

\oS J V ^ N 

~ASAp = (JTj ASAT = kT 0 (H) , (125) 

< 126 > 

35ss=(g^msir=»3i(g / ( 127 ) 

ABA,1= (m) VtN {w) P ASAT =^ {m) p - (128) 

Table [l] summarizes all fluctuation relation derived for the NVT ensemble. 


9.5.2 NpT ensemble 


Now suppose that the volume of the system can also fluctuate while N =const. The 
stiffness matrix is now 2x2 and its Kn component is still given by Eq. (11121b For the 
remaining components we have 


/ dS\ fdp_\ 

= J_(dT\ 1 V^/T.iV WJ P 

112 kT 0 \dVj SN kT 0 fdS_\ Nkc v ’ 

V J v,n 

where we used the Maxwell relation 

/ dS \ / dp \ 

V ) T,N V ) V.N 


and 

l/(9p\ 1 / \ 2 1 f dp \ 

22 = ~kT 0 {dVj StN = lSTkc~ v \&f)~ W 0 \dVj TN ’ 

where we used the identity 

/ dp 5 / dp \ T / dp \ 2 

\dh/ SAf di / 7 v It:, &I J VN 


(129) 


(130) 


(131) 


( 132 ) 
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(A E) 2 = NkTyCy 


(AS) 2 = Nkc v 
(AE)(AS) = NkT 0 c v 


(AT) 2 = 


kT[ 

Nc v 


(A^) 2 


(A p) 2 


kTl(ds ) 2 
Nc v \dT J p 

kT[(dp \ 2 

Nc v \dTJ p 


ATAp 


ApAT 


kT[fdp\ 
Nc v \dTJ p 

kT[ (dp\ 
Nc v \dT J p 


ApAp 


kTl(dp\ fdp\ 
Nc v \dT) p \dTj p 


AS AT = kT 0 


AEAT = kT* 


ASAp = kT 0 



p 


ASAH=kT 0 (^^ 

=kTS (If) f 


Table 1: Complete set of fluctuation relations for a simple fluid in the NVT ensemble. 






















Note that the derivatives (dp/3 V)t,n and ( dp/dV)s,N are related to the isothermal and 
adiabatic moduli, respectively. 

Thus, the stiffness matrix of the fluid takes the form 


/ 


K = 


Nkc v 


with the compliance matrix 


dp 

OT 


\ 


dp 

dT 


dp \ 2 Nc v / dp 

fff)~^\dVj TN J 


(133) 


K-^-kTof^- 

\ dp 


( f dp\ 2 Nc v f dp 
{dT/ \dV 


T,N 


T,N 


V 


dp 

TFT 


dp\ \ 
OT 




cm 


We are ready to calculate the fluctuations. To find the entropy and volume fluctuations, 
we apply Ea. ffT(j3l) : 


(AST = = NkC - Wig) g )’= Nkc r , 


dV 


(Ah) 2 = = - kT ; 


dP / T,N 

dV\ (dp\ lrT1 fdV\ 


ASAV =- kT ° I wL(w)r kT ° Uv 


(135) 

(136) 

(137) 


where c p is the specific heat at constant pressure. Fluctuation relations involving energy are 
calculated by Eq. dllOp : 


= Nktfc, - KI* g 


S,N 


T,N 


dp \ Po 


dT/ p T 0 


(138) 


AEAS = 


dE\ 


as J v , N Kl ' + {dv) S:N k ™ 


2 ( dV\ f dp 


NkT 0 c v - kT/ ( — ) ( ^ 


\ dp J 


. dT 

T,N \ u± / p 


dp 

dT 


Po 

T 0 


(139) 
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AEAV = 


dE\ 


I KE 1 + ( ^77 I K. 


( dE\ 


L 12 


) V,N 

(SL 


= - kT { 3 2 | AA ] 


\dV J 

dp 
dT 


L 22 


S,N 


Po 

To 


( 140 ) 


Turning to fluctuations of intensive parameters, we use Eqs. (11061) and (11071) to obtain 

(141) 


(AT) 2 = (kTo) 2 K u = 


kT[ 

Nc ’ 


(Ap) 2 = {kToYK 22 = -kTo 


dp 

dV 


T 0 (dp 


T,N 


Nc„ V dT 


kT( (dp 


ATA P= -(kT oY K 12 = J £[£) . 


(142) 


(143) 


To derive fluctuation relations involving the chemical potential, we treat it as a function of 
T and p. Using Eq. dllip . 


(Ap,) 2 =(kT 0 ( 


dp 

dT 


K n + 


A( k 22 - 2 (A) (A 


dp J 


T 


dT J \ dp 


K 


12 


(kT 0 ) 2 [s 2 0 K u + VqK 22 + 250 ^ 0 -^ 12 ] 


kT[ 

Nc v 


Vo 


dp 

dT 


2 


where 


so = - 

is the equilibrium entropy per particle and 


so 


dp 

dT 


kT °< (w 


T,N 


••=( 1 ), 

is the equilibrium volume per particle (u 0 = 1/po)- Similarly, 


(144) 


(145) 


(146) 


A pAT = (kToY 


dp\ 

dT), 


K u - 


dp 

dp 


K 


12 


kT[ 

Nc v 


Vo 


dp 

dT 


— s 0 


(147) 


ApAp 


(kToY 


A\ k 12 + (A 


dT ) 


dp 


K. 


-kT 0 v 0 I 


kT[(dp_ 
dVj TN Nc v \dT 


22 


VO 


dp 

dT 


So 


(148) 

(149) 
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Equation (I105j) gives the “cross-fluctuations' 


AS AT = -AVAp = kT 0 , 


ASAp = 0, 


(150) 

(151) 

(152) 


AVAT =0. 

The remaining “cross-fluctuations” are computed in a manner similar to Eqs. fllOSD - flllljh 


ASAp = 


dp 

dT 


AS AT = —/cT 0 s 0 , 


ABAT =( I 


AS AT = kT), 


V,N 


AEAp = I gy J AVAp = kT 0 po, 


(153) 

(154) 

(155) 


AE ^ = 1 If) 


V,N 


dp\ 


dE\ f dp 


dV J 


S,N 


dp 


AS AT + — AVAp = kT 0 (p 0 v 0 -T 0 s 0 ), 


dp 


AVAp = —— AVAp = — kT 0 v 0 . 

\dp J T 

The fluctuation relations for the NpT ensemble are summarized in Table [2l 


(156) 

(157) 


9.5.3 pVT ensemble 

We now turn to the grand-canonical ensemble, in which the volume of the system is fixed 
while the energy and the number of particles can fluctuate. The equilibrium temperature 
T 0 and the chemical potential p 0 are fixed by the reservoir. The independent fluctuating 
variables are S and N, with the conjugate thermodynamic forces P\ = T, and P 2 = p. 
Energy is a dependent extensive parameter. The chemical potential and pressure are also 
dependent parameters treated as functions of T and p. 

All fluctuation relations can be obtained from those for the NpT ensemble by simply 
swapping the variables V N, p —$■ —p. In particular, the derivative ( dp/dT) p is replaced 
by —( dp/dT) p and ( dp/dV)r,N by — V~ 1 (dp/dp)x■ Furthermore, it easy to verify that s 0 
is replaced by sq/vq and vq by 1/vq. The stiffness matrix becomes 


K 


N 0 kc v 


V©, 


(dpS \ \ 

\dTj p 

dp \ 2 Po c v fdp\ 
dT) p T 0 \dp) T ) 


(158) 
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( A s)> = Nkc-mi%) $ 


_ fdV\ 

ASAV = kT 0 l^—j 


p,N 




dV\ 


A/xA p = -kT 0 v 0 ( J-p 


T,N 


kri [dp 


Nc v V dT 


vo 


AS AT = kTn 


AVAp = -kT 0 


dp 

dT 


(A E) 2 = NkT$c v 
\ dp / t n 


dp \ po 


dT) p T 0 


ASAp = 0 


AEAS = NkT 0 c v 

MwUf), 


dp 

dT 


Po 

T 0 


AEAV = ~ kT ° 2 {%) T , N 


' dp \ _ Po_ 
KdTj 0 To 


_ UT 2 

(A T) 2 = -A 

NCy 


■kT 0 


(A pY 

_ UT i 

ata ^n)_ 


' dp \ 
\dV ) T)JV 


Tp {dp N 
Ac, V<9T, 


<? fdp\ 
.dT ) p 




ut 2 


^0 


<9p 

AT 


s o 


AV AT = 0 


ASAp = —kT 0 so 


AEAT = kT^ 


AEAp = kT 0 p 0 


AEAp = kT 0 (p 0 v 0 - T 0 s 0 ) 


AVAp = —kT 0 v o 


Table 2: Complete set of fluctuation relations for a simple fluid in the NpT ensemble. 
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(AS) 2 = 


ASAN 


(AN) 2 

(A Ey = 

+kTS ( 


AEAS 



AEAN 


(AT) 2 = 


(A/i) 2 = 


ATA/i 


(A p) 2 = 


ApAT 


Tabic 3: 


Pokc v V + kT 0 





ApAp 


kTp 

p 0 v 0 V 


kT{ / dp 
PoCvVqV \dT 




p 


AS AT = kT 0 


AN Ap = kTp 


PokTyCyV 



V 


p 0 kT 0 c v V 




Po 

n 


V 


_ kT 2 I ^ 
° [dp 


-l 


'(d_p\ 

po 

W) P 

'To 


V 


kTl 

P 0 CyV 


kT 0 

PoV 


Po 


dp 

dp 


+ 


dp 

dT 


kT 2 


dp' 

p 0 c v V \dT, p 


kTp {dp' + 


kT 2 


VpV \ dp J rp PpCyVpV 


s 0 + 77^ 


dp 

dT 


ASAp = 0 


ANAT =0 


ASAp = kTp— 
v 0 

AEAT = kT 2 


AEAp = kT 0 p 0 

_ ppp 

AEAp = - (pp + TpSp) 

v 0 

ANAp = — 

vo 


kT 2 s fdp\ 

poCyVpV So+ \dTj p 


Complete set of fluctuation relations for a simple fluid in the pVT ensemble. 


33 

















































with the inverse 


K~ l = kTaV 


dp\ 

dp) 


-i 


( (dp_ 

l dT 


V 


+ 


PoCy 

dp 

dT 


dp' 
dp / t 


dp\ \ 
dT 




(159) 


The results of the calculations are summarized in Table [3j 

The following scaling properties of fluctuations are evident from Tables [Q [2] and [3j 


• Covariances of extensive properties scale in proportion to the system size (N or V). 

• Covariances of intensive properties scale as inverse of the system size (IV -1 or V -1 ). 

• “Cross-fluctuations” between extensive and intensive properties are independent of the 
system size. 


The first of these properties explains why we could neglect the quadratic terms in the 
expansion for the reservoir entropy, justifying the reservoir approximation (j66j) discussed in 
Sec. 19.11 


10 Finite-reservoir ensembles 


10.1 General considerations 

We now revisit the system coupled to a reservoir discussed in Sec. 19.11 (Fig. |3j). The 
system and the reservoir are described by extensive parameters Xi,...,X n and X[,X^, 
respectively, with fixed total amounts = Xi + Xf. Suppose m of these parameters 
are allowed to flow back and forth between the system and the reservoir, whereas the re¬ 
maining (n — m) parameters are fixed. Both the system and the reservoir are assumed 
to be single-phase substances with different fundamental equations S = S(X 1 ,..., X n ) and 
S r = S r (X [,..., V”), respectively. The probability distribution of the fluctuating parameters 
of the system is given by Eq. (l64l) . which is repeated here for convenience: 


W(X u ...,X m ) = W m exp 




x exp 


SriX! - X u ...,X m - X m ) - S r (X, - ...,X m - X°J 

k 


(160) 


In Sec. 19.11 the entropy of the reservoir was approximated by the linear expansion (166|) 
at the equilibrium parameter set X°,..., X^. This linear approximation was justihed by the 
large size of the reservoir in comparison with the system. We will now lift, this assumption 
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and adopt a more general expansion, Eq. (l65j) . which includes quadratic terms: 


m 

S r (X 1 -X l ,...,X m -X m )-S r (X 1 -X° 1 ,...,X m -X° m ) = -^if(X,-X°) 

1=1 



-~k V A),(V - Xf)(X, - A'»), 

(161) 


i,j =1 


where 

1 ( d 2 S r \ 

(162) 


A ■ - 

k \dX[dXr) o 


is the stability matrix of the reservoir. As usual, the derivatives are taken at the equilibrium 
parameter set X°, ...,X^. 

Before analyzing fluctuations, let us examine the micro-state probability distribu¬ 
tion P(Xi ,..., X m ). Recall that the latter is obtained from W(X \,..., X m ) by dropping 
S(Xi, ...,X m ) in the first line of Eq. (11601) [see Eq. (TTCTT) in Sec. 19.2] . Thus, 


P = W m exp 
= A exp 


-S(X ?, ...,X°) + S r {Ad - X ls ...,X m - X m ) - S r {Ad - Xl ...,X m - X°) 


k 


m 1 m 

% £ F ° x > - ? £ A t< A * - x °)( x i - x 


i =1 


i,j =1 


(163) 


where A is a constant that can be found from the probability normalization condition. 
Thus, in addition to the usual linear terms FfX t appearing in the generalized canonical 
distribution (J72]h we now have a quadratic form representing the effect of the finite reservoir 
on the micro-state probability of the system. 

The matrix elements A£- scale as 1/N r . When the size of the reservoir increases to 
infinity, the quadratic terms in Eq. (11631) vanish and we return to the generalized canonical 
distribution (172]) . In the opposite limit when the reservoir becomes extremely small in 
comparison with our system, the probability distribution is dominated by the quadratic 
terms. The distribution becomes a very sharp Gaussian peak. Accordingly, all micro¬ 
states whose Xi are even slightly deviated from X? become extremely improbable. The 
system behaves as if isolated and Eq. (11631) approaches the micro-canonical distribution. In 
other words, Eq. (11631) interpolates between the generalized canonical and micro-canonical 
distributions; we can smoothly transition from one distribution to the other by varying the 
reservoir size and thus scaling the reservoir stability matrix (11621) . 

We now return to the parameter distribution function W(X i,...,X m ) and expand the 
entropy of our system in a manner similar to Eq. (1161 j) : 


S(Xi ,..., X m ) — S(X 


o 

!)•••) 



iil hl 

y f?(X, - X°) - -k A«(V - Xf)(Xi - .¥»), (164) 

i=l i,j =1 


with the usual stability matrix 























Inserting the expansions (1 161 j) and (I164j> in Eq. fllGOp . we obtain the multi-variable Gaussian 
distribution 


W (Ad, X m ) = W m exp 


A«(V - X°)(Xj - Xf) 


*,j=l 


(166) 


with the stability matrix 

A y = Ay + Ay. (167) 

We thus arrive at the remarkable result that the finite size of the reservoir does not 
change the Gaussian form of the parameter distribution of the system. However, the reser¬ 
voir does contribute to the matrix elements of the Gaussian distribution in an additive 
manner. We can, therefore, easily compute the mean-square fluctuations and covariances of 
the parameters Xi, ..., X m using the known Eqs. (18TT) and (1821) with A t] replaced by the full 
stability matrix A if 

AXiAXj = Aj- 1 , ij = 1,..., m, (168) 

(WW? = A^, i = (169) 


A number of other interesting relations can be derived. It follows from 


(11611) that 



FJ 



KAXj-x*), 

3 = 1 
m 

F° + kY,KA x i- x °)- 

3 = 1 


Eqs. (11641) and 


(170) 

(171) 


Averaging these equations over the distribution and using the fact that for a Gaussian 
distribution Xj = XJ, we obtain 

F[=T i = F°. (172) 


In other words, not only the most probable but also average values of the thermodynamic 
forces are the same in the system and in the reservoir. Furthermore, subtracting Eq. (11701) 
from Eq. (I171|) we have 

m 

FT-F^k^hiXi-X?). (173) 

i=i 

Starting from this equation, we can generate a set of fluctuation relations involving the 
difference ( F[ — F,). For example, multiplying Eq. (II731) by A Xj and averaging over the 
parameter distribution, we have 


{F[ — Fi)AXj = kSij, i,j — 1 , ..., m, (174) 

where we used Eq. (11681) . Likewise, multiplying Eq. (1173ft by (FJ — Fj) and averaging over 
the distribution, we obtain 

(FT - Fj)(FJ - Fj) = k 2 Aij, i,j = 1, ...,m, (175) 
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where we used Eq. (ll74p . In particular, 


(F[ — Fj) 2 — k 2 An, i — 1, ..., m. (176) 

These equations describe fluctuations of the differences between the values of the intensive 
parameters inside our system and in the reservoir. 

We can also compute fluctuations inside our system. To this end, we rewrite Eg. (117011 
as 

m 

AF l = -kJ2K AX P- (I 77 ) 

p =i 

Multiplying this equation by A Xj and averaging over the distribution, we obtain the co- 
variances 

m 

AFiAXj = -kYA ip A-?, i, j = 1,..., m. (178) 

p= i 

Similarly, multiplying Eq. (1177ft by A Fj and averaging over the distribution, we have 

m 

AFiAFj = k 2 Y ApA^A” 1 , = ( 179 ) 

p,q= i 

m 

(A Fi) 2 = k 2 Y A ipKq~Kii i = ( 18 °) 

p,q= i 

As expected, in the limit of an infinitely large reservoir when A tJ becomes identical to A,j, 
Eqs. (|178p . (1179ft and (I180p reduce to Eqs. (l85P (186P and (1871) for the generalized canonical 
ensemble. 


10.2 Application to a simple fluid 

To demonstrate the finite-ensemble formalism, consider a fixed volume V of a simple 
fluid embedded in another (generally, different) simple fluid serving as a reservoir. This 
situation was already discussed in Sec. 19.5.31 in the context of the grand-canonical ensemble. 
This time, however, the reservoir will be treated as an infinite source of heat (thermostat) 
but a finite source of particles. We thus have two fluctuating extensive properties, — E 
and X 2 = A, with the conjugate thermodynamic forces Fi = 1/T and F 2 = —p/T. 

Because the reservoir has an infinite heat capacity, we have Aj] = 0. Assuming for 
simplicity that the chemical potential of the reservoir /i r is temperature-independent, we 
additionally have A r 12 = A 2 i = 0- Thus, the reservoir stability matrix is 


A r . = 
o 


0 0 
0 A 


22 


(181) 


with a single nonzero parameter A r 22 . 

The micro-state probability distribution for our system can be found from Eq. (1163ft : 


P = A exp 


E AoN _1 2 

kT 0 2 22 v 


(182) 
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As a side note, a similar equation can be obtained for a closed system coupled to a thermostat 
with a finite heat capacity. In this case, the role of the number of particles is played by the 
energy, so that the additional term contributed by the thermostat is quadratic in (E — E 0 ): 


P = A exp 


E 

kTn 


-A r u (E-E„f 


(183) 


This situation is known in the literature as the Gaussian ensemble 47], [50]. By analogy, 
formula (I182j) can also be referred to as the Gaussian ensemble distribution. 

Equation (I182j) can be rewritten the form 


P = A exp 


E — pN 

kTn 


where 


with 


p = a — bN 


a — Po + A^ATqAo, 


b = -A r 22 kT 0 . 


(184) 

(185) 

(186) 
(187) 


Equation (11841) looks similar to the grand-canonical distribution, except that the fixed chem¬ 
ical potential of the reservoir p 0 is replaced an “effective” chemical potential p linearly depen¬ 
dent on the current number of particles N. In the limit of A 22 —> 0 (infinitely large reservoir), 
this dependence disappears and we return to the standard grand-canonical distribution with 
the chemical potential /i 0 . 

Equations (11841) - (11871) can be implemented in Monte Carlo computer simulations us¬ 
ing an algorithm similar to that of the grand-canonical ensemble [lj but with a chemical 
potential adjustable “on the fly”. For binary mixtures, this approach was implemented in 


the semi-grand canonical mode as “variance-constrained” Monte Carlo [54, |55|. A simi¬ 


lar semi-grand canonical Monte Carlo method, called feedback Monte Carlo, was recently 
employed for interface free energy calculations by the capillary fluctuation approach [2l| . 
Related Monte Carlo algorithms were earlier developed for closed systems coupled to a finite- 
capacity thermostat. The initially proposed “dynamical ensemble” method [59] developed 
for such systems later evolved into more sophisticated and computationally more efficient 
iterative Monte Carlo schemes based on Eq. (j!83j) [Hit , 53’. 

We now return to the analysis of fluctuations. To simplify the calculations, we will neglect 
the temperature dependence of the chemical potential, treating it solely as a function of the 
particle density p. Following the calculation methods presented in Sec. 19.51 it is easy to 
obtain the stability matrix of the fluid in the entropy scheme, 


/ 1 


A ij — 


p 0 c v kT^V 


~p0 Po c vTq 


~p0 

dp 

dp 


+ pI 


(188) 
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and its inverse 


Ar. 1 = kT 0 V 


-l 


/ 


d p 

dp y rp 


Accordingly, 


/ 1 


Kij p 0 c v kT 2 V 


ip--- 

\ 



Ho 


(189) 

Ho 

1 ) 



~H0 


\ 



dp 


-p 0 p 0 c v T 0 ( —) + pi + p 0 c v kT 2 A r 22 V 


(190) 


and 


~ Kij fdp 


kTnV 


— J + kToV A 22 
dp) T 


( 3 \ 

—J + pi + poc v kT) A 22 V p 0 


Ho 


\ 

!/ 


(191) 


We can now compute the mean-square fluctuations of the energy and the number of 
particles: 


(A E) 2 = A- = 


Po c v kT)V ( ^ ) + kT 0 V (p 2 0 + p 0 c v kT 2 A r 22 V) 


dp 

dp 


(AN) 2 = A 2 ) = 


+ kT 0 V A 22 
t 

kT 0 V 


dp 

dp 


(192) 


(193) 


+ kT 0 V A 2 


22 


Similarly, for the covariance AEAN we have 


AEAN = At) = —— 
' dp 


kT 0 p 0 V 


(194) 


—j T + kn v A ' 22 


Fluctuations of the chemical potential could be found directly from Eq. (118011 . Instead, 
we will take a shortcut by taking advantage of the fact that p = p(p) and therefore 


A p = 


dp\ AN 

d^) T ~V~' 


(195) 


Taking a square of this equation and averaging over the distribution, we obtain 

kT 0 [ dp 


dp\ (AN) 2 _ V \dp J T 


^=(lf P J T v . 


dp 

dp 


(196) 


+ kT 0 VA 


22 
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where we used Eq. (ll93p . Applying the same method, 


A ^ N = [f 


dp 

dp 


dp\ (AN) 2 




dp\ 

dp J T 


ApAE = 


V 


AEAN 

V 


tF) +kT 0 VA r o 


dp 

dp 


HokTo ( 'A 


(197) 


22 


dp 

dp 


(198) 


— I + kT 0 VA 2 , 


22 


Continuing along this line, all other covariances can be readily calculated. 

Consider limiting cases of the above relations. When A r 22 —> 0 (the reservoir is an infinite 
source of particles), we exactly recover all results for the grand-canonical ensemble listed 
in Table [3] When A 22 —y oo (the reservoir cannot supply or absorb particles), the system 
behaves as virtually closed (NVT). Accordingly, the energy fluctuation (1192)1 reduces to 


(A E) 2 = P okTjc v V 


(199) 


in agreement with the NVT result (Tabled]). All other mean-square fluctuations and co- 
variances appearing in Eqs. (11931) - (11981) tend to zero as 1/A 22 . As evident from Tabled! in 
the NVT ensemble they are proportional to (dp/dT) p and thus generally not zero. In the 
present calculations, they vanished because we neglected the temperature dependence of p 
for the sake of simplicity. More accurate but rather tedious calculations taking into account 
the temperature dependence of p exactly recover the fluctuation relations from Table CD when 
A 22 —>■ oo. 

It is interesting to note that at the limit of normal thermodynamic stability when 
(i dp/dp)x —> 0, the canonical fluctuation (AN) 2 diverges to infinity (Table [3]). However, 
when the system is coupled to a finite-size reservoir with a positive coefficient A 22 > 0, the 
latter stabilizes the fluid and its compositional fluctuations remain finite and equal to 1/A 22 , 
see Eq. (ll93p . The finite size of the reservoir smooths the singularity. 


11 Fluctuations of grain boundary properties 


GBs are interfaces separating regions of the same crystalline solid phase (called grains) 
with different crystallographic orientations. GBs can strongly impact physical and me¬ 
chanical properties of crystalline materials, especially in nano-structured systems. While 
thermodynamics and kinetics of GBs has been studied extensively for several decades 0, 


it is only recently that fluctuations of GB properties have become the subject of dedicated 
research, primarily by atomistic computer simulations 01221- 


In preparation for the discussion of GB fluctuations, we will give a briefly account of GB 
thermodynamics focusing on a plane GB in a binary substitutional solid solution. 
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Figure 5: Two grains with different crystallographic orientations separated by a plane grain 
boundary (GB). The shaded areas bounded by dashed lines mark a region containing the 
GB and a comparison region inside one of the grains containing the same total number of 
atoms. These two regions are selected for calculations of GB excess quantities. 


11.1 Grain boundary thermodynamics 


Consider a bicrystal (i.e., a system of two grains) enclosed in a rectangular box with 
rigid impermeable walls (Fig. [5j). Since the system is closed, the total number of particles 
N is fixed. We will assume that the grains possess appropriate crystallographic symmetries 
such that the mechanical stresses existing in the bicrystal do not create a difference in strain 
energy densities in the grans that would cause GB migration. 

Before discussing GB thermodynamics, we will first specify thermodynamic properties 
of the homogeneous single-crystalline solid solution forming the grains. We postulate the 
following fundamental equation of this solution in the energy representation^ 

E g = E a (S„V a ,N‘,Ni). (200) 


Here, the suffix g indicates that the properties refer to a region in one of the grains, such 
as the shaded layer shown in Fig. [5j TVf is the number of atoms of species 2 and N 9 is the 
number of lattice sites in the region. Vacancies are neglected, so that the number of atoms 
of species 1 is Nf = Njj — iVf. The differential of Eq. (12001) is 


dE g = TdSg 


pdV g + pdN 9 + MdN 9 , (201) 

where T = dE g /dS g is temperature, p = —dE g /dV g is pressure (negative of the stress 



Eg = TEg - P Vg + ^ + MiVf . 


( 202 ) 


15 


Generally, the fundamental equation of a stressed solid depends on the deformation gradient relative 


to a chosen reference state [2J, |61|, |62|. In the particular case considered here, the cross-section of the 


solid is fixed and only deformation in the normal direction is permitted. Under such conditions, only the 
normal component of the deformation gradient is relevant and varies as a linear function of the volume. It 
is convenient to represent this component by the volume itself. 
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Combining Eqs. (1201ft and (I202p we obtain the Gibbs-Duhem equation 


- SgdT + Vgdp - N°d<p - N 9 2 dM = 0. 


( 203 ) 


To develop GB thermodynamics, we define the GB excess of any extensive property X 
by 

X = X - X g , (204) 

where X is the value of the property for a bicrystalline layer containing the GB and X g is the 
value of the same property for a homogeneous grain layer containing the same total number 
of atoms. Examples of such layers are shown by shaded regions in Fig. [5] The choice of these 
layers is totally arbitrary and does not affect X as long as their boundaries are unaffected 
by the GB or the walls. According to this definition of excess, a bicrystalline layer can be 
thought of as composed of two subsystems: (1) the reference grain layer containing the same 
total number of atoms and (2) the excess system representing the GB and described by the 
excess properties S, E, V, etc. 

Just as the grain thermodynamics has been formulated above starting from a postulated 
fundamental equation (12001) . we will postulate a fundamental equation of the GB in the 
form 

E = E(S, V, N 2 , A). (205) 

Note that Afi is not an independent variable because 

Ni + N 2 = 0 (206) 

by our definition of excesses. The excess energy (1205jl is a homogeneous function of first 
degree in all four arguments. Applying the Euler theorem (66] we have 


^ dE ~ dE~ dE ~ dE A 
E = S + — V + —^ N 2 + —A 


dS dV dN 2 

On the other hand, the differential of Eq. (12051) is 


8A 


(207) 


dE ~ dE ~ dE ~ dE 
dE = —dS + — dV + ~^rdN 2 + — dA 


dS 


dV 


dN 2 


dA 


(208) 


Combining Eqs. (1207ft and (1208ft we obtain the GB version of the Gibbs-Duhem equation: 


r( dE ~ dE ~ , OE A dE n 
Sd—^r + V d-^ + N 2 d—^~ + Ad^-r = 0. 


dS 


dV 


dN 2 


dA 


(209) 


Next, we find the conditions of thermodynamic equilibrium between the GB and the 
grains. To this end, we consider a variation of the total energy (E + E g ) of a bicrystalline 
layer at ax fixed composition, volume and entropy. In equilibrium, this variation (dE + dE g ) 
must be zero [24|. Here, dE g is given by Eq. (120lft and dE by Eq. (1208ft . The variation is 
subject to the following constraints: 


dN x + dN{ = 0, 


( 210 ) 
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dN 2 + d/Vf = 0, (211) 

dv + dVg = 0 , ( 212 ) 

dS + dS g = 0 , (213) 

as well as dA = 0. Furthermore, by adding Eqs. (|210j) and (12111) and taking into account 
Eq. (12061) . it follows that the number of sites in the reference grain region is conserved: 

dNf = dNf + dN% = 0. Inserting these constraints into the total energy variation, we arrive 

at the equilibrium condition 


dE \ ~ dE \ ~ dE 

—— — T ) dS + [ ———|- p ) dV + [ — — 
dS / \dV \ dN, 


M]dN 2 = 0. 


(214) 


Because the variations dS, dV and dN 2 are independent and can be positive or negative, 
the coefficients multiplying them must be zero. We thus obtain the equilibrium values of 
the derivatives of E: 

dE . 

— T 0 , 


= M 0 , 


= -Po- 


As usual, the subscript 0 labels the equilibrium state. 
The remaining derivative 

( dE 

70 " (si 



(215) 

(216) 
(217) 


(218) 


is called the equilibrium GB free energy. The precise thermodynamic meaning of this prop¬ 
erty will be discussed later. 


11.2 Quasi-equilibrium grain boundary states 

So far, we have only specified the derivatives of E with respect to S, N 2 , V and A 
for equilibrium states of the system. Namely, such derivatives are equal to the equilibrium 
temperature, diffusion potential, negative of pressure in the grains and the equilibrium 
GB free energy, respectively. For the analysis of GB fluctuations, we must consider these 
derivatives without assuming equilibrium between the GB and the grains. Such derivatives 
form a new set of variables T, p, M and 7 defined as follows: 


dE 

~d§ 




dE 

aw 


= M, 


(219) 

( 220 ) 
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( 221 ) 


B E 
BV 
BE 
dA 


= ~P, 


= 7- 


( 222 ) 


These new variables can be interpreted as, respectively, the local temperature, negative 
of local pressure, local diffusion potential and the GB free energy when the GB is not in 
equilibrium with the grains. Indeed, by postulating the fundamental equation (12051) we 
implicitly assumed that the GB follows this equation even in the absence of equilibrium 
with the grains. One can think of such GB states as being disconnected from the actual 
grains and equilibrated with imaginary grains with a different set of intensive parameters 
T, M and p. In terms of the variables (I219j) - (|222|) . the differential form of the fundamental 
equation (I205p becomes 

dE = TdS - pdV + MdN 2 + ydA. (223) 


The existence of quasi-equilibrium states in which the GB obeys a fundamental equation 
without being in equilibrium with the grains is an important supposition of the interface 
fluctuation theory. As discussed in Secs. [4] and [5j the entire thermodynamic theory of 
equilibrium fluctuations is built on the assumption that all fluctuated states are states of 
quasi-equilibrium. The theory discussed here is an application of this principle to GBs. 

Suppose the bicrystal has been initially equilibrated and then the GB deviates from this 
initial equilibrium to a nearby quasi-equilibrium state where it is no longer in equilibrium 
with the grains. Denoting the equilibrium excess entropy, volume and number of atoms So, 
Vo and iVi), respectively, the Gibss-Duhem equation (I209j) for this variations becomes 


~ dE ~ BE ~ 0 , BE A BE 
Sod —- + Vn d —— + N~ d —-—b Ad—— 


BS 


BV 


bn 2 


BA 


= SodT - V 0 dp + N%dM + Ady = 0. 
This equation can be rewritten in the form 

dy = --A-dT + -ydp --p-dM, 

A A A 


(224) 


(225) 


which will be used later. 

In the particular case when the states of the grains also vary in such a manner that the 
entire system undergoes an equilibrium process, Eq. (l225p becomes 

dy 0 = -^dT 0 + |dp 0 - ^dM 0 . (226) 

This is one of the forms of the Gibbs adsorption equation for a fixed GB area. This equi¬ 
librium equation should not be confused with Eq. (l225p describing fluctuations in which the 
GB deviates from equilibrium with the grains. 

We can now understand the meaning of the parameter y. Let us rewrite Eq. (l207p in the 
form 

E = TS — pV + MN 2 + y A, (227) 
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from which 


(228) 


7 = 


E-TS + pV- MN 2 
A ' 


As was shown by Gibbs 24j, an expression of the form E — TS+pV — YjiPiNi is the reversible 
work required for the creation of a new system with extensive properties E, S, V and IV* 
by drawing the energy and matter from an infinite reservoir with a temperature T, pressure 
p and chemical potentials /i,. Thus, the numerator in Eq. (12281) can be interpreted as the 
reversible work of GB formation from an imaginary infinite reservoir with the intensive 
parameters T, p and M. In the present case, the term p\N\ + p 2 N 2 reduces to MN 2 due 
to the constraint imposed by Eq. (12061) . In short, 7 is the reversible work of GB formation 
per unit area. Accordingly, the equilibrium GB free energy 70 is the reversible work of GB 
formation per unit area from a reservoir with the same intensive parameters To, po and Mq 
as in the actual grains: 

E-T 0 S+p 0 V-M 0 N 2 

7o --^-• (229) 

In other words, this is the reversible work of GB formation between two infinitely large 
grains. 


11.3 Grain boundary fluctuations 

We are now in a position to analyze GB fluctuations. We will compute generalized 
canonical fluctuations of GB properties treating the grains as an infinite reservoir cou¬ 
pled to the GB. Although the foregoing discussion utilized the energy scheme, the latter 
is impractical for GB fluctuations^ Instead, we will now switch to the entropy scheme in 
which the independent additive invariants fluctuating between the GB and the grains are 
X\ — E, X 2 — V and X 3 = N 2 . The conjugate thermodynamic forces are readily found 
from Eq. (l223j) : F\ = 1/T, F 2 = p/T and F 3 = —M/T. We can also rewrite Eq. ()225[) in 
terms of the new variables: 


dy = ^2 


1=1 


with the coefficients 


B\ = 


T 0 (E 0 — 7 o A) 


B 2 = 

b 3 = 


A 
T 0 Vo 
A ’ 

t 3 n 2 


A 


(230) 

(231) 

(232) 

(233) 


16 The energy scheme involves fluctuations of the excess entropy S , which is not accessible by experiments 
or atomistic simulations. In the entropy scheme, S is replaced by the more accessible excess energy E. 
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In the Gaussian approximation to fluctuations, the distribution of the excess parameters 
X t is given by Eq. (TT9]i with the stability matrix 


Aij k 


( ( d 2 S 
dE -2 

d 2 S \ 
dEdV) 

( d 2 S 


\\8EdN 2 


d 2 S 

dEdV 

'd 2 s' 


dV 2 

d 2 S 
dVdN o 


d 2 S 
dEdN 2 

d 2 S 

\dVdN 2 

d 2 s' 

dN 2 


(- 


\ 


/ 


(234) 


The elements of this matrix represent equilibrium excess properties of the GB. For example, 


d 2 S 

])W 2 



(dE (t))„ T 0 Mc„’ 


(235) 


where 

i (oe\ 

c °-a\w) o 

is the GB heat capacity (more accurately, the excess heat capacity per unit GB area com¬ 
puted at fixed V/A and N 2 /A). Similarly, 



d 2 S 

dV 2 


d_ 

dV 


Po 


ToVokr T 2 V 0 a 


where 



is the isothermal compressibility of the GB and 


o 


1 




o 


(237) 


(238) 


(239) 


is the GB thermal expansion coefficient. 

Equilibrium GB properties can be computed directly by atomistic simulations. For 
example, Rt can be found by computing the excess GB volume as a function of pressure 
at a fixed temperature followed by numerical differentiation |67|]. However, analysis of 
GB fluctuations offers a more efficient approach in which the entire set of equilibrium GB 
properties can be extracted from a single simulation run using fluctuation relations. Indeed, 
the excess properties E, V and N 2 are readily accessible by atomistic simulations. By 
computing their covariances, the inverse stability matrix can be found from the fluctuation 
relation 
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AX,A Xj = A"- 1 , i,j = 1,2, 3. (240) 

The latter can be then inverted to obtain the stability matrix A tJ and thus the equilibrium 
GB properties such as c v , Rt, <5, etc. Furthermore, knowing the stability matrix, fluctuations 
of thermodynamic forces can be predicted from the fluctuation relations 

AFiAFj = k 2 Aij, i,j — 1,...,3. (241) 


Knowing these covariances, straightforward algebraic manipulations can be applied to find 
(AM) 2 , (AT) 2 and (A p) 2 , as well as the covariances AM AT, AMAp and ApAT. As a 
general rule, 



where the transformation matrix dPi/dF m is 

(-T 2 p 0 T 0 -M 0 T 0 \ 



V 0 0 -T 0 / 


(242) 


(243) 


To our knowledge, the proposed calculation scheme has not been implemented. 

As an interesting application, one can predict the equilibrium mean-square fluctuation 
(Ay) 2 of the GB free energy, where Ay = y — y 0 . The fluctuation form of Eq. (l230p is 


3 

Ay = J2 BiAF i- 

i =1 


(244) 


Taking a square of this equation and averaging over the distribution, we obtain 

3 

XF=E BiBj AFiAFj. (245) 

i,j =1 

Applying Eq. (12411) for AFiAFj, this equation finally becomes 

3 

(Ay) 2 = k 2 Y, B iBjAij. (246) 

i,j =1 

This calculation requires the knowledge of A i3 and the equilibrium values of the excess GB 
energy, volume and segregation as ingredients for the parameters 5,. The GB free energy 
is usually treated as a static property. An estimate of its fluctuation (Ay) 2 would provide 
a critical assessment of this approximation. 
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Figure 6: Sharp interface model of GB pre-melting. The pre-melted GB is modeled by a 
liquid layer of a width w bounded by two sharp interfaces endowed with thermodynamic 
properties of actual solid-liquid interfaces between bulk phases. 

12 Fluctuations in pre-melted grain boundaries 

In this section we address another fluctuation topic associated with GBs. At high temper¬ 
atures approaching the melting point of the solid, many GBs become atomically disordered 
and can develop a liquid-like structure reminiscent of a thin liquid him j(30|. Our goal is 
to describe equilibrium thermodynamic properties of this him and derive the distribution 
function of its width. The analysis is based on the sharp interface model of pre-melting 
(Fig. [6]). The grains are formed by the same binary solid solution phase as discussed in 
Sec. CU 

12.1 Thermodynamic background 

12.1.1 Equilibrium state of pre-melted grain boundary 

We will consider the same bicrystal in a rigid box as discussed in Sec. [HI but this 
time we will adopt a particular model of GB structure, namely, a liquid layer of a certain 
width w bounded by two solid-liquid interfaces (Fig. [6]). We start by describing equilibrium 
thermodynamic properties of this boundary. 

As in Sec. [HJ thermodynamic properties of the solid solution forming the grains are 
fully described by the fundamental equation (I200p in the energy representation. The binary 
liquid solution (index L ) is described by its own fundamental equation 

E l = E l {S l ,V l ,N*‘,NZ) (247) 


with the differential form 

dEi = T l cISl — PlcIVl T p\dN^ + /qd-ZV^''. (248) 

Here, T L and Pl are the temperature and pressure of the liquid and pi and /i 2 are the 
chemical potentials of the components. 

To describe thermodynamics of the solid-liquid interfaces bounding the liquid layer, we 
will first take a detour and consider a single plane interface between the bulk solid and liquid 
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phases. We can choose an arbitrary geometric dividing surface partitioning our system in 
two parts with volumes V g and Vl- We then compute the excess of any extensive property 
Wby 

X = X - X g - X L , (249) 

where X is the value of the property for the solid-liquid system with a volume (V g + Vl) 
and X g and X L are properties of the phases computed as if the volumes V g and Vj j remained 
homogeneous all the way to the dividing surface j24|. By this definition, V = 0. Because 
the densities of all properties are generally different on either side of the interface, the 
excess defined by Eq. (I249p depends on the exact position of the dividing surface within the 
interface region. We will adopt the rule by which we always place the dividing surface in 
such a manner that the excess of the total number of atoms N vanishes: N = 0. The 
fundamental equation of the solid-liquid interface is formulated in the form 

E = E(S,N 2 ,A), (250) 

with the differential form 

E) p 1 Ej 77 1 E\ rp 

dE = —-dS -I- —dN 2 + -r—rdA. (251) 

OS dN 2 dA K ’ 

Returning to the pre-melted GB modeled by a thin layer of the liquid phase (Fig. [6]), we 
will assume that the solid-liquid interfaces bounding this layer follow the same fundamental 
equation (1250ft as if they were interfaces between two bulk phases. This is, of course, an 
approximations. In reality, the liquid layer width is on the order of a nanometer and its 
thermodynamic properties in this narrow confinement are different from properties of bulk 
liquid. In addition, the solid-liquid interfaces separated by such a short distance can interact 
with each other due to their hnite thickness comparable with the width of the layer. The 
difference between the actual GB structure and this highly idealized liquid-layer structure 
is accounted for by introducing an interaction between the two solid-liquid interfaces called 
the “disjoining” interaction. The disjoining interaction (index d) is included in the model by 
adding an extra energy term Ed and an associated entropy Sd, and postulating the following 
fundamental equation of the disjoining interaction: 

E d = E d (S d ,V L ). (252) 

The dependence on Vl has been added to make the disjoining interaction a function of the 
distance between the solid-liquid interfaces, which for a fixed cross-section considered here, 
is proportional to Vl■ Both Ed and Sd are proportional to the GB area and have the same 
dimensions as the excess quantities E and S, respectively. Both Ed and Sd tend to zero if 
the width of the liquid layer increases to infinity (bulk solid-liquid system). 

To summarize the model, the GB is represented by a homogeneous liquid layer separated 
from the grains by two geometrically sharp solid-liquid interfaces to which we assign ther¬ 
modynamic properties of actual solid-liquid interfaces between bulk phases (Fig. [6j). These 
sharp interfaces coincide with the dividing surfaces separated by a distance w and interact 
with each other by disjoining forces that follow the fundamental equation (1252(1 . 

We will next derive the conditions of thermodynamic equilibrium by considering a vari¬ 
ation of state of a bicrystalline layer indicated in Fig. [6] by a shaded region with dashed 
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boundaries. During this variation, the layer remains isolated from the rest of the system 
and keeps a constant value of its entropy. Under these constraints, the variation of the total 
energy 


E — Eg + El + E + Ed 


(253) 


must be zero 24]. The terms of this equation are represented by the fundamental equations 


(I200p . (1247 ]) . (12501) and (1252(1 . respectively. Accordingly, 


dE = TdS g - pgdVg + (pdN% + MdN$ 

+ T^dSi, — PlEVl + pidN^ + p 2 dN 2 


dE ~ dE ~ dE,. 

T ~ dS + ~ dN 2 + ~ . dA 

8S dN 2 dA 

dE ddc , dE dj ^ n 

+ oc dSd + „ dV l — 0 
dS d dV L 

(254) 

subject to the following constraints: 


dS g + dSL + dS + dSd = 0, 

(255) 

dVg + dVi, = 0, 

(256) 

dN? + dNf + dN x = 0, 

(257) 

dN? + dN 2 L + dN 2 = 0, 

(258) 

dN\ T dN 2 = 0, 

(259) 

dN? + diN? - dN f = 0, 

(260) 


as well as cL4 = 0. Here, — p g is the component of the mechanical stress in the grains normal 
to the GB. 

Out of 14 differentials appearing in the right-hand sides of the Eqs. (12541) - (12601) . 7 are 
eliminated by the constraints, leaving 7 independent variations: 


+ f^r -PL+Pg^J dV L + ((p-n^dNf 

( 8E \ 

+ hi — P 2 j dN 2 = 0. (261) 

Equating the differential coefficients to zero, we obtain the following equilibrium conditions: 


T 




P°L = P° q + 


dE 0 

dVr 


(262) 

(263) 
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Mo — /J-2 A*i j 


( 264 ) 


dE \ 

m) o 

ifo = n\, ( 265 ) 

where To is the uniform temperature throughout the equilibrium system. As usual, the 
superscript 0 marks the equilibrium state. 

Equation (1263ft is the mechanical equilibrium condition, which can be written in the form 


P°L = P°c, + Pl 


where p d is the equilibrium value of the disjoining pressure defined by 


Pd = 


dE d 

dV L 


( 266 ) 


( 267 ) 


Thus, the liquid phase within the GB layer is subject to an additional pressure p d on top of 
the pressure p g existing in the grains. 


12.1.2 Calculation of the disjoining pressure 

Calculation of the equilibrium disjoining pressure p d can be simplified by reformulation 
of the equilibrium conditions in the density form. The fundamental equation (1200ft of the 
solid solution can be rewritten in the density form 

£g = £g( S g, Vg, Cg) , (268) 

where £ g , s g and v g are, respectively, the energy, entropy and volume per lattice site (which 
is same as per atom since we neglect vacancies), and c g is the site (atomic) fraction of species 
2. Accordingly, T = d£ g /ds g , p g = —d£ g /dv g and M = d£ g /dc g . From Eq. (12021) we have 

ip = £g-TSg+ PgVg ~ MCg, ( 269 ) 

while the Gibbs-Duhem equation (I203p becomes 

SgdT — Vgdpg + dp + CgdM = 0. (270) 

It is more convenient to express all properties of the grain as functions of the variable 
set ( T,p g ,c g ). The thermodynamic potential appropriate for these variables is 

^(T, Pg, Cg) =<f+ MCg = £g~ TgSg+PgVg. (27 1 ) 

Its differential form is 

du = —SgdT + Vgdpg + Mdcg , (272) 

which is easily derivable from Eq. (1270ft . Similarly, for the liquid solution we will use the 
variable set (Tl,Pl,cl) and, accordingly, the Gibbs free energy 

g(T L , p L , c l ) = £ L - T L s L + p L v L = (1 - c L )pi + c L p 2 (273) 
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with the differential 


(274) 


dg = -s L dT L + v L dp L + (/i 2 - Pi)dc L , 

where all densities are counted per atom. 

We now reformulate the equilibrium conditions (I264p and (I265p in the form 


M ( T o,P°n,c°) = H2 (T 0 ,p£,c£) - /ii(T 0 ,^,4), 


0 „0' 


„0 „0 


(275) 


W (T„, p% c°) = g(T 0 ,pi cl) + M(T„, p° f , - cl), (276) 

respectively. This is a system of two non-linear equations with respect to five equilibrium 
parameters (T 0 ,p° g ,p 0 L ,c g ,c 0 L ). It follows that our solid-liquid system has three degrees of 
freedom. Let us choose the independent variable set (T 0 ,p°,c°). For any given state of the 
grains defined by these parameters, we can solve Eqs. (I275p and (12761) to find the composition 
c° L and pressure p° L in the liquid layer, and thus the disjoining pressure p d = p° L — p° g . 

For practical applications, we will replace this exact calculation scheme by a simple 
approximation. At given temperature T 0 and pressure p° g , let c* and c* L be the solid and 
liquid compositions corresponding to the bulk solid-liquid coexistence (p d = 0). These two 
compositions must satisfy the equation^ 


M ( T o,p° g ,c * g ) = p 2 (T 0: p° g ,c * L ) -/ii (T 0 ,p° g ,c* L ), 
w ( T o ,P° q ,c* q ) = g{T 0 ,p° c* L ) + M(T 0 ,p°, c*)(c* — c* L ). 


(277) 

(278) 


Note that the pressure continuity condition p° L = p° g valid for a two-phase bulk system with 
a plane interface eliminates one degree of freedom, leaving our system with two degrees of 
freedom in compliance with the Gibbs phase rule. For a fixed pressure Eqs. (12771) and 
(12781) define the solidus and liquidus lines, c*(T 0 ) and c* L (T 0 ), on the phase diagram of the 
system. 

We seek an approximate form of Eqs. (12751) and (12761) in the vicinity of the solidus line. 
As a measure of proximity of the grains to the solidus line at given T 0 and p g , we choose 
the diffusion potential deviation 


AM = M(T 0 ,p°, c°) — M(T 0} p 


c 


(279) 


from its solidus value M(T 0 ,p®,c g ). All functions appearing in Eqs. (l275p and (I276D can be 
then approximated by linear expansions in the small parameters AM, (c° g — c*), ( c° L — c * L ) 
and ( p° L — p° g ). The linearized form of these equations becomes 


AM 


f 3(^2 ~ 

V dc L A 



( d(p 2 - Mi) 
V dp L 


Pd, 


* 


(280) 


doj 

dc„ 


(4 - c * 9 ) 


dg 

dc r , 


LT - C*r) + 


dg 

dp L 


P° d 


+ M(T 0 ,p°, c!)(cj - c° L -c* + c* L ) + (c! - c^)AM, (281) 


17 We will not go into details, but Eas. (1277|l and (12781) can be given a geometric interpretation of a common 
tangent construction, M being the slope of the common tangent to the curves ui(To,p g ,c) and g(To,p®,c) 
as functions of composition c at fixed To and p° g . 
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where the asterisk indicates that the derivatives are taken in the bulk-coexistence state 
(AM = 0). All zero-order terms have canceled out by Eqs. (|277p and (12781) . The equations 
obtained are simplified by utilizing the relations 


fdu_ 

\9cg 


' dg N 

) =M(T 0 ,pl,cl), 

(282) 

dg > 
dp L J 

) =4, 

* 

(283) 

PiV 

- j 

\ —2 ^ 

) = V L ~ V L , 

* 

(284) 


dp L 

where v l L = dni/dpi, are partial molar volumes of the components in the liquid phase, with 
the asterisk indicating that these volumes refer to the bulk-coexistence state. Using these 
relations, Eqs. (l280p and (128ip become, respectively, 


AM = 


8 -% ^ ') (4 - 4) + («T - ®I*)pS, 


o = 44 + (4 - 4)A m. 

Equation (I286D immediately gives us the equilibrium disjoining pressure, 


(285) 

(286) 


n Cr C n 

P° d = j ^am, 


(287) 


while Eq. (l285p can be solved for the deviation of the liquid-layer composition from the 
liquidus composition c* L : 


, 2 * 


1* 


Ct = 


A M - (vl* - v l M 
~ Pi) 
dcr, 


1 - 


V L - V L 


{cl ~ C* ) 


d(P2 - Pi) 
dCr 


-AM. 


(288) 


In many systems, the difference between the partial molar volumes of the components is 
small in comparison with vl and the latter equation can be simplified to 


Ct - Ct = 


d(n 2 - Pi) 
dc T , 


-AM. 


(289) 


For a single-component system, AM loses its significance. The appropriate measure 
of proximity of the system to bulk melting is then the difference (To — T m ) between the 
temperature of the bicrystal and the bulk melting point T m at a given pressure p° g . The 
equilibrium conditions (I275p and (I276p reduce to one equation, 


u{To,P° g ,0) = g{T o ,p° L ,0) } 


(290) 
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which for bulk equilibrium becomes 


o) = g(T m ,p° o). 


(291) 


Linearizing Eq. (l290p in the small parameters (T 0 — T m ) and p d we have 


(292) 


from which 


pS = - T„) = 2L_(T 0 - T m ), 

Ut -L mV T 


J L - L mv L 

where H m = ( s* L — s*)T m is the enthalpy of melting per atom 


(293) 


12.2 Analysis of fluctuations 

To describe fluctuations in the pre-melted GB, we treat the grains as an infinite reservoir 
and apply the canonical fluctuation relation (j98lh The system coupled to the reservoir is 
the GB possessing the energy 

Eqb = El + E + Ed. (294) 

We will initially consider the full set of fluctuating parameters, which will later be reduced 
to a smaller set. Namely, we will initially treat all arguments of the fundamental equations 
(1247]h (I250p and (I252p . except for the fixed GB area, as the fluctuating parameters. In the 
energy scheme, the fluctuating parameter set is 


Z = (S L} S,S d ,V L ,N^N^N 2 ). 

The thermodynamic forces conjugate to these parameters are 

p_ /rp dE 3E d SE d dE 

i L 'ds' as d ’ PL+ dv L ’ t ‘ u,J ' 2 ’afij' 

The reversible work of GB formation (1991) takes the form 

7 


Tl=E GB -E° aB -Y,P°(Z-Zf), 


(295) 


(296) 


(297) 


i=l 


where the thermodynamic forces are taken in the state of equilibrium. Using the equilibrium 
conditions f!262p - f|264p . we obtain 

n = (E l - E° l ) + {E - Eo) + (E d - E° d ) 

-T 0 (S l - S° L ) - US - So) - T 0 (S d - S° d ) 
+pl{V L -Vl)-p° d {V L -Vl) 

—piUi - U°) - p°(iY 2 L - N™) - M 0 (N 2 - iV 2 °). (298) 
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Next, we will make some approximations that will reduce the number of fluctuating 
parameters. The difference ( El — E^j) appearing in Eq. (l298|) will be approximated by its 
linear expansion around equilibrium: 

El-E° l = T„(S l - Si) - pKV l - V°) + pKNf - AT") + p°(N^ - iV 2 "). (299) 

Likewise, (E — E $) will also be evaluated in the linear approximation: 

E-E 0 = T 0 (S-S 0 ) + M 0 (N 2 -N°). (300) 

The terms related to the disjoining interaction will remain intact. Inserting the linear 
expansions (12991) and (I300p in Eq. (12981) . several terms cancel out and we are left with 

n = {E d - E°) - T 0 (S d - S °) - p° d (V L - V L °). (301) 

Although Eq. (l301l) contains variations of three variables, only one fluctuating parameter 
is independent. To see this, we apply the Legendre transformation with respect to the 
disjoining entropy to obtain the thermodynamic potential 

T(T, V L ) = E d (S d , V L ) - TS d . (302) 

This thermodynamic potential is called the disjoining potential and its differential form is 

dT = —S d dT + p d dV L . (303) 

The latter equation suggests another definition of the disjoining pressure: 

Pd= {w) T (304) 

Thus, Eq. (l30ip can be rewritten in the form 


n = T(T 0 , V L ) - T(T 0 , V° L ) - p° d (V L - V®), (305) 

showing that, at a given temperature T 0 , 1Z is solely a function of V dj . In other words, V dj is 
the only independent fluctuating parameter. By neglecting higher-order terms in Eqs. (l299p 
and (13001) . we have suppressed fluctuations of all other independent parameters. 

The probability distribution of the liquid volume Vl is given by Eq. (|98l) . which becomes 


W{V L ) = W m exp 


TZ 


kTn 


= W m exp 


T(r 0 ,14) - T(r 0 , vl) - p^v L - V g) 

kT, 


(306) 


Remembering that the cross-section of the bicrystal is fixed, the finally arrive at the distri¬ 
bution function of the liquid-layer width w: 


W(w ) = W m exp 


-A 


y?(T 0 , w) - ^{Tp, wp) - p° d (w - wp) 
kTn 


(307) 
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where Wq = V^/A is the equilibrium width of the layer and ^ = 4//A is the disjoining 
potential per unit area. Note that this distribution depends on the GB area A: the larger 
the area, the sharper is the peak of the distribution around w 0 . The disjoining pressure can 
be rewritten as pd = (d^ / 8w)t■ Its equilibrium value p d appearing in Eq. (|307|) can be found 
from the bulk thermodynamic properties using Eq. (l287j) for a binary alloy and Eq. (l293|) for 
a single-component system. Knowing p d and the distribution of re at a given temperature 
T 0 , one can invert Eq. (13071) and extract the disjoining potential ip(T 0 ,w). 

The GB width distribution (13071) in conjunction with Eqs. (l287p and (I293p for the equi¬ 
librium disjoining pressure answers the goal of this section. For single-component systems, 
an equation similar to (13071) was derived in previous work [16 - 18|, [68| and utilized for calcu¬ 
lations of disjoining interactions in Ni GBs by atomistic simulations 16-118]. In the present 
paper, this distribution has been extended to binary systems. 

As a further approximation, we can expand A (To, w ) around the equilibrium width keep¬ 
ing only linear and quadratic terms: 


fp(T 0 ,w) = ip(T 0 ,w 0 ) + p° d (w - w 0 ) + i ) O-^o) 2 , 


(308) 


where the second derivative is taken at equilibrium. This results in the Gaussian distribution 


W(w ) = HA exp 


A 


8 2 ^ 


2kTr\ \dw 2 


(w 


w 0 f 


L o / 0 

from which the mean-square fluctuation of the GB width is 


(w - w 0 ) 2 = 


kT ( 


0 


A(d 2 il>/dw 2 )o 


(309) 


(310) 


This estimate shows that the root-mean-square fluctuation of the liquid-layer width scales 
with the GB area as l/\/A 

According to the foregoing derivation, the GB width w appearing in the sharp-interface 
model of GB premelting (Fig. [6]) has the meaning of the distance between two dividing 
surfaces of the solid-liquid interfaces. Such dividing surfaces are characterized by a zero 
excess of the total number of particles [V = 0, A^ = 0, see excess definition (12491) ]. It 
should be, therefore, understood that the alternate definitions of w adopted in atomistic 
simulations 16-[l8] constitute approximations to the strict definition. 


13 Concluding remarks 

The thermodynamic theory of fluctuations presented here expands the postulational ba¬ 
sis of the traditional classical thermodynamics by incorporating equilibrium fluctuations. 
The main additional element of the expanded theory is the fluctuation postulate expressed 
by Eq. (lT5]) [or equivalently, Eq. (fT9l) ]. Due to its statistical nature, this postulate signifies a 
break with the determinism of classical thermodynamics. The statistical component intro¬ 
duced by this postulate propagates through the entire logical structure of thermodynamics. 
The equilibrium state of a system is no longer a state but a distribution over states. Ther¬ 
modynamic properties become distributions aroimd the average values which are measured 
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by macroscopic experiments. Similar expansions of classical thermodynamics were previ 


ously discussed by Tisza and Quay [39| and Callen |26|, [271, |38|] and can be traced back to 


Einstein 40]- The fluctuation postulate formulated in this paper describes fluctuations in an 
isolated system and stems from the microcanonical distribution in statistical mechanics. In 
this sense, our approach is different from that of Tisza and Quay [39] and Callen 26, 27, 38 
and is more aligned with Einstein’s proposal (4o| . 

The fluctuation postulate relies on the concept of non-equilibrium entropy. Without 
a proper definition of non-equilibrium entropy, we cannot talk about its increase during 
equilibration or behavior during fluctuations. In some sense, non-equilibrium entropy is 
defined already in classical thermodynamics. The thermodynamic equilibrium conditions of 
classical thermodynamics are formulated through the entropy maximum principle { 24 J, 0- 


28]. In this principle, the actual equilibrium state of the system is compared with possible 


virtual states in which the system is imagined to be partitioned into isolated equilibrium 
subsystems; their entropies are then summed up and compared with the equilibrium entropy 
of the actual system. The expanded thermodynamics goes one step further: it assumes 
that such states are actually implemented during spontaneous fluctuations away from the 
equilibrium state. The term “quasi-equilibrium” introduced in this paper captures three 
defining properties of such fluctuated states: (1) non-equilibrium, (2) real (as opposed to 
virtual states of classical thermodynamics), and (3) composed of subsystems that can be 
treated as isolated (on a certain time scale) and each described by a fundamental equation. 
The quasi-equilibrium states are characterized by a set of internal parameters Aj. The 
fluctuation postulate specifies the distribution function of these parameters for an isolated 
system. 

The second law of thermodynamics is re-formulated to include fluctuations. In the refined 
formulation, it is only the entropy S averaged over the fluctuations that monotonically in¬ 
creases and reaches a maximum at equilibrium. On shorter time scales, the non-equilibrium 
entropy S incessantly fluctuates around its average value and can never accede a certain 
maximum value S (Figs. [2] and [3]). Thus, the average entropy S obtained by macroscopic 
measurements (performed on the time scale bro) is always somewhat smaller than S. 

In classical thermodynamics, a central role is played by the fundamental equation which 
encapsulates all equilibrium thermodynamic properties of the system. This equation, first 
discovered by Gibbs [0], looks exactly like our Eq. (J4]) but has a subtly different meaning. 
In the classical fundamental equation, the entropy S' is a static property defined through 
the maximum-entropy principle (e.g., Callen’s Postulate II [0. 0| or Tisza’s Postulate Phi 


28]) and measured experimentally. In the expanded theory, S appearing in Eq. (J4|) has 


the statistical-mechanical meaning of the maximum micro-canonical entropy of an isolated 
system. As noted in the previous paragraph, this entropy is slightly larger than the measured 
entropy S. In other words, the measured entropy S does not exactly satisfy the fundamental 
equation ((4]) of the present theory. This fact, of course, does not create any problems. 
It only reflects the possibility of different definitions of the fundamental equation for a 
fluctuating system. Since both S and S are well-defined for any given set of conserved 
parameters X\, ...,X n , S(X \,..., X n ) could also be taken as the fundamental equation [26I.127I. 
3 38j- We would then have to re-define the temperature, pressure and all other intensities, 
and to modify many other equations of the theory. In the end, taking the fundamental 
equation S(X \,..., X n ) as the point of departure would not lead to a logically consistent 
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thermodynamic fluctuation theory consistent with statistical mechanics P*1 The advantage 
our Eq.(j5J) with the entropy defined as above is that it leads to a simple and complete logical 
structure fully consistent with statistical mechanics. In the limit of an infinite large system, 
the fluctuations vanish and S and S become numerically equal. Accordingly, Eq.flT]) becomes 
the fundamental equation in the sense of classical thermodynamics. 

Similar considerations apply to canonical fluctuations. If the average values X\,X m 
of the fluctuating extensive parameters Xi,X rn are formally inserted in the fundamental 
equation (J4]) , we do not obtain the canonical entropy S. However, in the limit of a large 
system, all canonical averages do satisfy ([4]) and the latter becomes the classical fundamental 
equation. 

For the sake of simplicity, we have only considered a particular class of thermodynamic 
systems for which all arguments of the fundamental equation are additive invariants. In 
the future, the theory can be generalized to include pseudo-thermodynamic and quasi¬ 
thermodynamic variables 28]|. Pseudo-thermodynamic variables are neither additive nor 


conserved; an example is furnished by long-range order parameters in atomically ordered 
compounds. Quasi-thermodynamic variables can be additive but need not be conserved. 
They give rise to additional work terms in the total energy variation. However, the conjugate 
thermodynamic forces do not become spatially homogeneous when the system reaches equi¬ 
librium and thus cannot be considered as intensities. Examples include the total magnetic 
or electric-dipole moment of the system, the conjugate variables for which are, respectively, 
the magnetic and electric field. Elastic strains and stresses in inhomogeneously deformed 
crystalline solids belong to the same category. All such cases require a separate treatment. 

To facilitate practical applications of the theory, we have presented a regular procedure 
for calculations of equilibrium fluctuations of extensive parameters, intensive parameters 
and densities in systems with any number of parameters. The proposed formalism has 
been demonstrated by deriving the complete set of fluctuation relations for a simple fluid 
in three different ensembles. The results are summarized in Tables [0 El and [3] and may 
present a reference value. It should be noted that these calculations treat all fluctuations in 
the Gaussian approximation. Greene and Callen 26, 0 proposed a more general method 
that does not rely on the Gaussian approximation and enables calculations of higher-order 
correlation moments such as AX, : AXj AX fc , AX* AXj AX*. AX/, etc. Such moments are 
rarely important and were not discussed in this paper. 

The proposed fluctuation formalism has been applied to solve two problems related to 
GBs in binary solid solutions. First, we developed a set of equations relating the fluctua¬ 
tions of excess GB properties, readily accessible by atomistic simulations, to equilibrium GB 
characteristics such as the excess elastic modulus, excess heat capacity and others. These 
equations offer a computationally efficient approach to computing the entire set of such 
properties from a single simulation run. Furthermore, the covariances between the fluctua¬ 
tions generate a set of thermodynamic relations that can be tested by simulations to verify 
the fundamental understanding of GB thermodynamics. Implementation of this program is 
presently underway. 


18 As already mentioned, Greene and Callen [37| resolved these difficulties by making an additional ap¬ 
proximation that essentially identifies the canonical and micro-canonical entropies. This identification is 
not satisfying for a rigorous fluctuation theory. The distinction between the two entropies was discussed in 
Sec.liOl 
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Second, we have extended the sharp interface model of GB pre-melting to binary alloys 
and derived the distribution of the fluctuating liquid-layer width w. This distribution has re¬ 
cently been applied to calculate the disjoining potential for several GBs in Cu-Ag alloys {23|. 
The disjoining potential is a key GB property allowing predictions of different pre-melting 
scenaria, such as the continuous melting, thin-to-thick transitions or abrupt melting after 
over-heating |(39j]. In the version of the pre-melting model considered here, we essentially 
treat the grains, the liquid layer and the solid-liquid interfaces as parts of the same infinite 
reservoir. We only account for fluctuations of w due disjoining interaction. All thermody¬ 
namic relations are in place to develop a more general model that would include fluctuations 
of intensive properties of the liquid layer and/or the solid-liquid interface properties. This 
generalization could be the subject of future work. 

Finally, the approach demonstrated here for GBs can be readily extended to other inter¬ 
faces. It would be interesting to study the fluctuations of excess properties of solid-liquid 
and especially solid-solid interfaces to determine their excess elastic and thermal properties 
and test some of the basic equations of interface thermodynamics. 
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